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INTRODUCTION 


The  following  annual  progress  report  documents  work  at  the  University  of 
Minnesota  Heat  Transfer  Laboratory,  under  AFOSR  sponsorship,  on  the  topic  of 
heat  transfer  from  gas  turbine  airfoil,  end-wall  and  internal  passage  surfaces. 
The  research  is  divided  into  the  following  subtopics:  curvature  effects, 
including  transition  and  film  cooling  and  airfoil  heat  transfer,  including 
portions  near  the  endwall.  Progress  on  each  is  discussed  separately. 


A.  The  Effect  of  Convex  Curvature  on  Heat  Transfer  and  Hydrodynamics  of 
Turbulent  Boundary  Layers  —  Including  the  Process  of  Recovery  from 
Curvature 

Three  experiments  were  conducted  to  investigate  the  effect  of  streamwise 
curvature  over  the  range  from  6/R  -  0.01  to  0.0*1  and  to  study  the  effect  of 
free-stream  turbulence  on  fully  turbulent  curved  boundary  layers.  An  attached 
fully  turbulent  boundary  layer,  grown  on  a  flat  plate  is  introduced  to  a 
convex  wall  of  constant  radius  of  curvature  followed  by  flat  recovery  wall. 


1 .  Objectives 

The  effect  of  varying  the  radius  of  curvature  on  curved,  two-dimensional 
boundary  layers,  including  the  recovery  process,  was  investigated  to  further 
understand  the  transport  mechanism  and  to  extend  the  data  base  for  prediction 
model  development.  Two  different  free-stream  turbulence  intensity  cases  were 
employed  to  study  the  effect  of  the  free-stream  turbulence  intensity  on  the 
curved  boundary  layer.  Significant  additions  to  the  previous  data  base  are: 

a.  The  recovery  process,  with  heat  transfer,  was  investigated  for  cases  of 
milder  curvature  than  in  previous  studies. 

b.  The  radius  of  curvature  was  varied  from  case  to  case  in  the  same  study. 

c.  Curvature  effects  for  cases  of  6/R  between  those  of  mild  and  strong 
curvature  were  investigated  in  the  same  facility. 

d.  Free-stream  turbulence  intensity  effects  on  the  curved  boundary  layer 
were  investigated  for  the  first  time. 


2.  Progress  and  Accomplishments 

The  important  observations  from  the  three  experiments  were  published  in 
You,  Simon  and  Kim  (1986a),  You,  Simon  and  Kim  (1986b)  and  Kim  and  Simon 
(1987).  Some  details  of  the  turbulent  heat  flux  measurements,  not  given  in  a 
previous  progress  report,  follow: 

The  effects  of  streamwise  convex  curvature,  recovery,  and  free-stream 
turublence  intensity  on  the  turbulent  transport  of  heat  and  momentum  has  been 
investigated.  A  special  three-wire  hot-wire  probe  based  on  the  work  of  Blair 
and  Bennett  (1984)  has  been  developed  for  this  purpose.  Two  cases  with 
free-stream  turbulence  levels  of  0.68J  and  2.0%,  taken  in  the  facility 
described  earlier  with  6/R  -0.03  (a  moderate  strength  of  curvature)  were 
compared.  A  summary  of  the  results  follows. 

1.  Profiles  of  u'v',  t* ,  u't*  and  v't’  are  reduced  by  curvature.  An 
asymptotic  profile  was  achieved  very  early.  Recovery  in  the  weak 
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portion  of  the  profile  oooure  quickly,  with  the  profiles  often 
overshooting  the  flat  plate  (upstream  of  the  curve)  values. 

2.  A  reversal  of  sign  in  the  shear  stress  is  seen  for  the  high  TI  case 
(TI-2.0?)  for  y/6  >  0.5.  This  behavior  was  observed  by  other 
researchers  in  strongly  curved  flow  with  lower  turbulence  intensity. 

3.  Beoause  the  measured  turbulent  Prandtl  number  profiles  show  considerable 
scatter,  one  can  only  conclude  that  the  effect  of  curvature  is  no  more 
than  about  20?.  But  the  data  appear  to  show  a  decrease  by  the  end  of 
recovery  from  curvature.  This  does  not  refute  Prt  values  deduced  in 
earlier  studies  from  mean  velocity  and  temperature  profiles,  but  does 
not  provide  support  for  the  20?  increase  noted  by  You,  et  al.  (1986a). 
Negative  turbulent  Prandtl  numbers  were  seen  in  the  wake  for  the  high  Tl 
case,  indicating  a  breakdown  in  Reynolds  analogy. 

A.  Curvature  effects  dominate  turbulence  intensity  effects  for  the  cases 
studied. 

The  reader  is  referred  to  Kim  and  Simon  (1987)  for  more  details. 

B.  The  Concave  Wall  Study 
1 .  Objectives 

A  major  weakness  in  the  design  of  gas  turbines  lies  in  the  ability  to 
predict  the  location  of  transition  from  laminar  to  turbulent  flow  on  the 
blades.  Graham  (1979)  stated  that  a  35?  error  in  the  prediction  of  heat 
transfer  coefficient  results  in  a  wall  temperature  error  of  56 °C,  and  an  order 
of  magnitude  decrement  in  blade  life.  Since  as  much  as  50-80?  of  the  blade 
can  be  covered  by  flow  undergoing  transition,  and  since  the  skin  friction  and 
convective  heat  transfer  increase  significantly  during  the  transition  process, 
the  ability  to  prediot  the  momentum  and  heat  transport  in  this  region  is 
crucial . 

Some  of  the  parameters  known  to  influence  boundary  layer  transition  are 
free-stream  turbulence,  acoustic  disturbances,  surface  vibration,  surface 
roughness,  streamwise  acceleration,  blowing  and  suction,  boundary  layer 
separation,  compressibility,  body  forces,  and  streamwise  curvature.  The 
effect  of  streamline  curvature,  pressure  gradient,  and  free-stream  turbulence 
are  investigated  in  the  present  study.  The  effect  of  the  other  variables  and 
various  other  uncontrollable  parameters  are  minimized  by  using  a  unique 
flexible  test  wall,  e.g,  Wang  (1984)  enabling  all  experiments  to  be  run  on  a 
single  rig. 

The  study  is  to  begin  with  a  baseline  case  on  a  straight  wall  with  no 
streamwise  acceleration  and  low  turbulence  intensity.  All  subsequent  runs 
will  be  compared  to  this  case  so  that  the  effects  of  the  above  three 
parameters  cm  transition  can  be  isolated.  The  proposed  experimental  runs  are 
summarized  below: 
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Case 

TI 

R(cm) 

Acc. 

Description 

1 

0.68 

• 

No 

Baseline  Case — Plane  surf..  Low  TI 

2 

0.68 

-180 

No 

45  deg-bend,  Low  TI,  No.  acc. 

3 

0.68 

-90 

No 

45  deg-bend,  Low  TI,  No.  acc. 

4 

2.0 

m 

No 

Plane  surface,  High  TI,  No  acc. 

5 

2.0 

-180 

No 

45  deg-bend,  High  TI,  No  acc. 

6 

2.0 

-90 

No 

90  deg-bend.  High  TI,  No  acc. 

7 

m 

Yes 

Plane  surf.,  streamwise  acc.,  Low  TI 

8 

-180 

Yes 

45  deg-bend,  Streamwise  acc.,  Low  TI 

9 

2.0 

m 

Yes 

rlane  surf.,  streamwise  acc..  High  TI 

The  above 

ac./dec. 

runs  will 

be  with 

constant  B  values  (self-similar). 

Special  oases  nay  be  added  where  the  acc. /dec.  values  simulate  the  blade 
pressure  side  and  the  turbulence  intensity  is  elevated  above  2%,  simulating 
blade  conditions. 

The  experimental  program  consists  of  measuring,  as  a  minimum,  the 
following  quantities. 

1.  Mean  and  fluctuating  components  of  streamwise  velocity.  This  is 
measured  using  a  pitot  tube  in  the  heated  flow,  and  by  a  horizontal 
hot-wire  in  isothermal  flows. 

2.  Mean  temperature  profiles.  The  thermocouple  probe  described  by  Wang 
(1984)  will  be  used. 

3.  Local  Stanton  number.  One  hundred  thirty  thermocouples  embedded  in  the 
test  wall  are  for  this  purpose. 

4.  Shear  stress  profiles  and  profiles  of  the  fluctuating  component  of 
cross-stream  velocity  in  isothermal  flows  where  the  boundary  layer  is 
sufficiently  thick.  A  cross-wire  probe  will  be  used. 

5.  Transverse  velocity  correlation  coefficient  in  the  free-stream  as  well 
as  profiles  of  the  autocorrelation  of  streamwise  veoocity.  Two 
horizontal  hot-wire  probes  will  be  used. 

If  the  boundary  layers  are  thick  enough,  measurements  of  the  three  components 
of  velocity  along  with  their  cross-correlations  will  be  made  with  a  triple 
wire  where  the  flow  displays  considerable  three-dimensionality.  Streamwise 
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vorticity  measurements  using  a  probe  similar  to  that  described  by  Wallace 
(1 984 )  are  also  being  considered.  _  _ 

A  probe  to  measure  the  apparent  turbulent  heat  flux  u't'  and  v't'  in 
two-dimensional  boundary  layers  similar  to  that  described  by  Blair  and  Bennet 
(1984)  has  been  developed.  The  reader  is  referred  to  Kim  and  Simon  (1987)  for 
details.  The  utility  of  this  probe  on  the  concave  wall  will  depend  on  the 
two-dimensionality  of  the  boundary  layer. 

If  there  is  a  formation  of  counter-rotating  longitudinal  vortices  on  the 
ooncave  wall  prior  to  transition,  the  flow  will  be  strongly  three-dimensional. 
It  is  therefore  necessary  to  know  where,  in  relation  to  the  vortices,  the 
measurements  are  taken.  Flow  visualization  is  thus  necessary. 

The  main  method  of  visualization  is  by  using  a  sheet  of  cholesteric 
liquid  crystal  bonded  to  the  surface  of  the  concave  wall.  The  crystals  change 
color  with  temperature,  enabling  the  temperature  field  to  be  mapped.  Also, 
with  a  constant  heat  flux  boundary  condition,  lines  of  constant  color 
correspond  to  lines  of  constant  heat  transfer  coefficient.  By  heating  the 
wall  slightly,  the  variation  in  wall  temperature  caused  by  the  longitudinal 
vortices  may  be  seen.  Lower  temperature  lines  correspond  to  downwash  between 
the  vorticies,  while  higher  temperature  lines  correspond  to  the  upwash.  This 
method  of  visualization  enables  one  to  monitor  the  formation  and  growth  of  the 
vortices  as  well  as  determine  their  spacing. 

Two  other  methods  of  visualization  are  being  considered.  Infrared 
scanning  of  the  concave  wall  would  also  enable  mapping  of  the  surface 
temperature  field.  An  IR  scanner  is  currently  available.  The  other  method  of 
visualization  being  considered  is  the  smoke-wire  method.  A  thin  wire  covered 
with  oil  is  stretched  across  the  test  section,  and  resistance  heating  of  the 
wire  vaporizes  the  oil,  producing  "smoke".  Visualization  of  the  vortices  is 
made  possible  due  to  lumping  of  the  smoke  at  the  upwash.  The  electronic 
circuit  which  drives  a  smoke-wire  system  has  been  built.  A  helium  bubble 
generator  is  also  available  for  visualization,  as  are  high  speed  movie 
cameras . 


2 .  Status 


The  test  wall  to  be  used  in  this  study  has  been  built.  A  schematic  of  it 
is  given  in  Fig.  1.  The  design  is  similar  to  that  described  by  Wang  (1984), 
with  the  exception  of  using  a  lexan/liquid  crystal  composite  instead  of  the 
stainless  steel  and  the  3M  P-19  film.  The  lexan  was  necessary  for  many 
reasons.  The  smoothness  of  the  lexan  eliminates  the  possibility  of  early 
transition  due  to  waviness.  It  also  permits  the  machining  of  a  leading  edge 
onto  the  wall.  A  measurement  program  in  which  the  emissivity  of  the  liquid 
crystal  was  measured  eliminated  the  need  for  the  P-19  reflective  film  that  had 
been  used  for  radiation  control.  The  thermal  conductivity  of  the  lexan/liquid 
crystal  was  also  measured. 

Installation  of  the  test  wall  has  been  completed.  Qualification  of  the 
test  wall  with  regards  to  transition  location  was  performed  by  heating  the 
wall  and  visualizing  transition  using  the  liquid  crystal.  Transition  was 
assumed  to  occur  at  the  location  where  the  liquid  crystal  first  changes  color. 
This  corresponds  to  the  hottest  wall  temperature  or,  since  the  wall  heat  flux 
is  uniform,  the  location  of  lowest  heat  transfer  coefficient.  Various 
parameters,  such  as  the  leading  edge  suction  flow  rate,  the  distance  from  the 
nozzle  exit,  eto.  were  optimized  such  that  transition  occurred  as  far 
downstream  as  possible  for  a  given  free-stream  velocity  and  turbulence 
intensity.  The  outer  flexible  wall  was  adjusted  such  that  there  was  no 
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pressure  gradient  along  the  wall.  The  Reynolds  number  based  on  displacement 
and  momentum  thickness  at  the  beginning  of  transition  were  measured  to  be  1920 
and  737,  respectively.  The  free-stream  turbulence  intensity  was  0.M3J.  A  plot 
of  Reg,  vs.  free-stream  turbulence  intensity  for  the  present  study  is  shown 
on  Fig.  2,  and  is  seen  to  be  in  excellent  agreement  with  the  data  of  previous 
researchers. 

Further  qualification  of  the  test  wall  will  be  performed  using  momentum 
and  energy  balances.  Programs  to  do  this  are  currently  being  written. 

C.  Modeling  of  Curved  Turbulent  and  Transitional  Flows 

Previously,  experiments  on  boundary  layer  transition  had  been  executed  in 
the  Heat  Transfer  Laboratory  (Wang,  Simon  and  Buddhavarapu  (1985)  and  Wang  and 
Simon  (1985)).  In  addition  to  this,  fully-turbulent  boundary  layer 
experiments  have  been  made  as  discussed  in  the  above  section  (You,  Simon  and 
Kim  (1986a),  You,  Simon  and  Kim  (1986b)  and  Kim  and  Simon  (1987)).  These 
sets,  all  taken  in  our  lab  with  the  same  experimental  techniques,  give  us  a 
unique  opportunity  to  evaluate  some  of  the  models  typically  employed  in  gas 
turbine  heat  load  evaluation;  in  particular,  empirical  transition  modeling  and 
empirical  correction  for  convex  curvature  effects  to  the  mixing  length 
hypothesis  turbulence  model.  The  results  of  this  work  are  given  in  Park  and 
Simon  (1987),  which  is  attached. 

1 .  Objectives 

a.  To  test  present  design  models  against  the  flat  plate  transitional  data 
and  the  convex-wall  curve-flow  data,  separately. 

b.  To  test  the  combined  model,  which  has  a  chosen  transition  model  and 
curvature  convection  to  the  turbulence  model,  against  the  curved-wall 
transition  data. 

c.  To  test  the  hypothesis  that  curvature  correction  need  by  applied  to  the 
transition  start,  length  and  intermittency  submodels. 

2.  Progress  and  Accomplishments 

The  work  is  documented  in  Park  and  Simon  (1987).  A  major  conclusion  is 
that  transition  models  can  be  found  which  "predict"  the  transition  data  but 
prediction  from  the  various  models  which  were  chosen  for  review  disagree 
somewhat  from  one  to  another  and  that  a  model  chosen  based  on  one  case  does 
not  "predict"  another  case  particularly  well.  A  most  useful  finding  was  that, 
when  free-stream  turbulence  intensity  (T.I.)  values  are  around  2t  (and  it  is 
presumed  the  same  could  be  said  for  higher  T.I.  values),  the  above  hypothesis 
is  true  and  curvature  corrections  need  not  be  incorporated  into  the  transition 
start,  length  and  intermittency  submodels.  When  the  T.I.  is  0.7t  (and  it  is 
presumed  to  be  true  for  values  below  1$)  the  hypothesis  fails.  This  lends 
some  support  to  present  design  methods  which  ignore  the  curvature  effect  on 
transition  location. 

3.  Expected  Activities  for  Next  Year 

Further  mixing  length  hypothesis  modeling  will  remain  "on  hold"  until  the 
concave  transitional  data  sets  have  been  completed.  Modeling  with  more 
sophisticated  and  general  turbulence  closure  models  will  continue  under  the 
direction  of  Professor  S.V.  Patankar. 


D. 


Film  Cooling  of  Curved  Surfaces 

1 .  Objectives 

The  purpose  of  this  study  is  to  secure  basic  understanding  of  the 
processes  involved  in  film  cooling  on  curved  surfaces.  More  specifically,  the 
behavior  of  flow  emerging  from  a  row  of  round  Jets  on  the  convex  and  concave 
walls  of  a  curved  duct  is  being  examined  with  particular  attention  to  how  well 
these  jets  protect  the  walls  of  the  duct  from  the  free-stream  conditions. 
Recently,  an  investigation  into  the  effects  of  injection  rate  and  curvature 
has  been  concluded.  A  brief  summary  of  the  findings  is  given  in  the  next 
section.  In  the  following  phase  of  this  study,  the  importance  of  injection 
angle  will  be  examined. 

2.  Progress  and  Accomplishments 

Film  cooling  through  a  row  of  holes  over  a  convex  and  a  concave  surface 
has  been  studied.  The  holes  were  inclined  at  35  degrees  to  the  main  flow  and 
were  spaced  three  diameters  center-to-center .  The  effects  of  injection  rate 
and  strength  of  curvature  on  performance  were  studied.  A  foreign  gas 
injection  technique,  employing  the  mass  transfer  analogy  to  heat  transfer,  was 
used. 

Large  variations  in  lateral  profiles  of  local  effectiveness  on  the  convex 
surface  were  found.  At  all  but  the  highest  blowing  rates,  the  local 
effectiveness  was  generally  high  at  points  in  line  with  an  injection  hole  and 
low  half-way  between  two  holes.  Lateral  profiles  on  the  concave  surface  were 
comparatively  flat.  This  was  attributed  to  concave  instabilities  which 
enhanced  mixing  there. 

Stronger  convex  curvature  was  found  to  enhance  lateral  average 
effectiveness  at  low  blowing  rates,  while  at  high  blowing  rates,  strength  of 
curvature  seemed  to  diminish  in  importance.  Also  at  high  blowing  rates,  the 
centrifual  force  of  the  Jet  caused  the  jet  to  lift  off  the  surface.  Momentum 
flux  ratio  was  determined  to  be  the  parameter  which  influenced  when  lift-off 
would  occur.  Further  increase  in  blowing  rate  or  momentum  flux  ratio  caused 
the  jets  to  merge  and  blockage  of  the  mainflow  from  the  wall  resulted.  An 
increase  in  effectiveness  was,  therefore,  encountered. 

For  film  cooling  of  concave  surfaces,  momentum  flux  ratio  had  relatively 
little  influence  on  the  performance  of  film  cooling  on  concave  surfaces.  This 
was  probably  because  the  normal  momentum  of  the  Jets  worked  to  degrade 
performance,  while  the  tangential  momentum  (centrifugal  force)  caused  an 
improvement.  With  the  present  injection  geometry,  the  effects  cancelled  each 
other  except  Just  downstream  of  injection,  where  high  normal  momentum  promoted 
lift-off.  After  twelve  diameters  from  injection,  blowing  rate,  rather  than 
momentum  flux  ratio,  was  the  dominant  parameter. 

On  the  concave  wall,  lateral  average  effectiveness  behaved  very  similarly 
to  analytical  correlations  for  slot  injection  on  a  flat  plate.  The  magnitude 
of  the  performance  was  consistently  lower  than  the  slot-flat  plate 
correlations,  but  the  trends  were  very  much  the  same.  Excessive  lateral 
mixing,  promoting  blockage  of  the  mainflow  from  the  wall  was  the  reason  for 
this  effect. 
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3.  Expected  Activities  for  Next  Year 

It  is  known  that  jet  centrifugal  force  and  normal  momentum  are  two 
important  parameters  determining  the  behavior  of  a  jet;  increased  Jet 
centrifugal  force  enhances  performance  on  the  concave  surface  but  degrades 
performance  on  the  convex  surface.  Jet  normal  momentum,  on  the  other  hand, 
degrades  performance  on  both  surfaces.  Both  injection  rate  and  injection 
angle  significantly  influence  Jet  normal  and  tangential  momentum. 

From  the  previous  phase,  understanding  of  the  effect  of  injection  rate 
has  been  obtained.  The  present  thrust  is  to  quantify  the  effect  of  injection 
angle. 

E.  Mass  Transfer  on  a  Turbine  Blade 

Better  understanding  of  the  characteristics  of  heat  transfer  on  a  turbine 
blade  will  improve  the  design  of  the  gas  turbine  engine.  Although  past 
researchers  reported  numerous  results,  these  data  mainly  were  measured  in  the 
two-dimensional  flow  region.  Also,  these  investigations  were  not  able  to 
show  detail  due  to  conduction  problems.  Because  the  ratio  of  span  to  chord 
length  is  not  large  enough  to  ignore  the  effect  of  the  endwall,  data  at  the 
two-dimensional  region  are  of  marginal  usefulness.  Goldstein  and  Karni  (198^) 
reported  that  mass  transfer  on  the  front  portion  of  the  cylinder  in  a  narrow 
span  next  to  the  endwall  is  increased  90  to  700  percent  over  the 
two-dimensional  region.  Further  investigation  of  heat  transfer  on  a  blade  in 
the  three-dimensional  region  is  necessary.  Flow  in  a  turbomachine  is  so 
complex  that  no  computational  calculation  cam  predict  it  accurately.  Direct 
measurement  therefore  becomes  essential  for  both  development  of  computational 
models  and  improvement  of  the  engine  design.  In  the  present  study,  a 
naphthalene  sublimation  technique  is  used.  Sublimation  rates  on  two  hundred 
and  fifty  locations  on  both  surfaces  of  a  blade  are  measured.  The  wind  tunnel 
and  a  linear  cascade  of  six  blades  used  in  the  study  are  shown  in  Fig.  3. 

1 .  Objectives 

First,  measurements  on  a  blade  in  the  two-dimensional  region  are 
conducted  to  compare  with  other  studies.  Because  the  Taylor-Gortler  vortex 
structure  might  occur  on  the  concave  surface  of  the  blade,  measurements 
searching  for  such  a  vortex  system  and  its  effects  on  heat  transfer  are  made. 
Finally,  measurements  will  be  conducted  near  the  end  wall  and  the  effect  of 
the  secondary  flow  system  on  heat  transfer  from  the  blade  surface  will  be 
carefully  investigated. 

2.  Progress  and  Accomplishments 

Presently,  the  data  acquisition  system  has  been  built  and  some 
preliminary  tests  have  been  made.  It  is  expected  that  most  of  the 
measurements  will  be  taken  at  the  end  of  this  summer. 

a)  Static  pressure  coefficients  are  measured  on  the  blade  surfaces 

and  the  approaching  velocity  is  measured  12.9  cm  upstream  of  the  blade. 
Values  of  pressure  coefficient  on  the  concave  side  and  convex  side  are  shown 
in  Fig.  M  and  Fig.  5,  respectively.  Values  of  X  are  measured  from  the 
stagnation  point  of  the  blade  and  S  is  the  chord  length.  Distribution  of 
pressure  coefficient  on  the  concave  side  is  quite  uniform.  The  nonuniformity 


of  the  convex  side  mainly  results  from  the  higher  pressure  at  X/S  -  0.43  which 
might  cause  the  suction  side  horseshoe  vortex  to  lift  away  from  the  wall. 

b)  Data  Acquisition  System 

The  flow  chart  of  the  data  acquisition  system  is  shown  in  Fig.  6.  The 
main  part  of  this  data  acquisition  system,  a  four-axis  table,  is  shown  in  Fig. 
7.  The  construction  of  this  table,  motor-controller,  and  interface  between 
all  individual  parts  have  been  in  progress  over  the  past  year.  The  software 
to  control  the  stepper  motors  and  LVDT  for  scanning  a  blade  surface  has  been 
tested.  Some  preliminary  tests  were  conducted,  i.e.  the  effect  of  the 
vibration  on  the  system  and  the  repeatability  of  measurements  on  a  solid 
cylinder. 

3.  Future  Activities 

We  plan  to  perform  measurements  on  a  blade  in  the  two-dimensional  region 
for  different  Reynolds  number  flows.  The  effect  of  the  Taylor-Gorlter  vortex 
structure"blade  heat  transfer  will  be  investigated  at  this  time.  Measurements 
of  mass  transfer  in  the  near-end-wall  region  will  then  be  conducted. 
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Fig.  1  —  Cross-section  of  Heated  Wall  (not  to  scale) 
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Fig.  3  Wind  tunnel  and  test  blades 
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Fravlous  studios  have  demonstrated  that  avan 
•lid  streamwlse  curvature  has  a  marked  offset  on 
turbulence  structure,  heat  transfer,  and  skin 
friction.  This  effect  ms  further  Investigated  In 
the  present  experiment  uhere  two  cases  of  dif¬ 
ferent  radii  of  convex  curvature  were  tested.  An 
attached  fully-turbulent  boundary  layer,  grown  on 
a  flat  plate  ms  Introduced  to  a  constant  radius 
of  curvature  convex  Mil  followed  fay  a  flat  reco¬ 
very  Mil.  The  strongth  of  curvature,  the  ratio 
ef  boundary  layer  thickness  to  radius  of  cur¬ 
vature,  6/R,  ms  0.013  and  0.03  for  the  two  cases, 
loth  cases  Mrs  taken  with  the  saae  test  Mil, 
bent  to  two  radii  of  curvature  and  the  saae  cur¬ 
vature  entry  conditions.  In  the  curve,  Stanton 
nuabers  and  skin  friction  Mrs  reduced  by  10%  and 
20%  below  flat  plate  values  for  6/R  ■  0.013  and 
(1.03,  respectively.  Recovery  ms  slow.  Turbulent 
Frandtl  nuabers  Mre  Increased  by  20-25%  returning 
rapidly  to  flat-Ml!  values  on  the  recovery  Mil. 
Shear  stress  profiles  responded  quickly  to  cur¬ 
vature,  reducing  to  about  6SX  of  the  upstream  flat 
plate  values  for  6/R-O . 03 . 


1.  INTRODUCTION 

A  better  understanding  of  transport  processes 
for  flow  over  stroaawlse  curved  surfaces  Is 
required  for  aaqy  practical  applications.  One 
representative  application  Is  In  gas  turbines 
irfiere  accurate  prediction  of  the  heat  transfer 
rates  Is  critical  to  a  successful  gas  turbine 
design.  Because  curvature  has  a  significant 
Influence  on  turbulent  boundary  layers,  high  per- 
fomence  turboMchlnery  design  aodels  oust 
appropriately  Include  Its  effects.  Presently,  the 
experimental  data  base  Is  not  sufficiently 
coaplete  to  develop  and  qualify  such  aodels.  Heat 
transfer  with  sustained  weak  curvature, 

•/R  •  0.01,  can  be  predicted  with  the  corrector 
proposed  by  Bradshaw  [1],  but  this  aodel  falls 
«dwn  curvature  bocoaes  substantially  greater  than 
*/R  »  0.01.  The  recovery  process  on  a  downstream 
flat  Mil  and  the  effect  of  curvature  on  turbulent 
Frandtl  number  are  presently  not  adequately 
•ode led,  even  under  M«k  curvature  conditions. 

The  objective  of  the  present  curved  boundary 
layer  study  Is  to  broaden  the  data  base  for  aodel 
developMnt.  For  this  study,  strength  of  cur¬ 
vature,  l/R,  ms  changed  by  changing  the  radius  of 
curvature,  R,  holding  the  entry  boundary  layer 
thickness  essentially  constant  The 


flexible  Mil  design  of  the  present  facility 
allows  such  a  study. 

Previous  Mork^ 

In  a  comprehensive  survey  of  literature  docu¬ 
menting  the  effects  of  streamline  curvature, 
Bradshaw  [1]  pointed  out  that  streaallne  curvature 
effects  on  shear  stress  and  heat  transfer  are 
significant  even  for  cases  of  6/R-0.01.  Some 
Important  findings  published  after  the  1968 
Bradshaw  survey  are  discussed  next.  In  1973,  So 
and  Mailer  [2],  showed  that  normal  and  shear 
stress  profiles  develop  self-slallar  shapes 
shortly  downstream  of  entry  to  a  strongly  curved 
section.  In  1979,  Slaon  and  Moffat  [3]  published 
heat  transfer  results  for  strong  convex  curvature 
and,  for  the  first  tlae,  a  flat  recovery  m11. 
Though  the  effect  of  curvature  ms  substantial , 
they  found  little  difference  between  cases  with 
6/R  ■  0.0S  and  0.08.  Recovery  of  Stanton  number 
differed  from  that  of  skin  friction  coefficient 
Indicating  a  curvature  offect  on  turbulent  Frandtl 
number.  Moyle,  et  al.  [4]  measured  local  heat 
transfer  coefficients  with  6/R  *  0.01;  120  boun¬ 
dary  layer  thicknesses  downstream  of  the  start  of 
curvature,  the  Stanton  number  ms  20%  less  than 
the  f lat-wal 1  correlation  on  the  convex  surface 
and  33%  greater  on  the  concave  surface. 

In  1982  and  1984,  Gibson  et  al.  [S,  6,  7] 
concluded  that,  with  6/R  *  0.01,  the  Stanton 
number  responds  wore  rapidly  than  the  skin  fric¬ 
tion  coefficient  to  curvature  change.  The 
Increase  of  turbulent  Prandtl  number  ms  21%.  The 
heat  transfer  results  of  Simon  and  Moffat  [8]  and 
the  fluid  aechanlc  measurements  of  61111s  and 
Johnston  [9]  showed  the  effects  of  strong  stream- 
wise  curvature  and  recovery  from  curvature  on  a 
downstream  flat  surface.  Shear  stress  profiles 
within  the  curved  region  for  two  strengths  of  cur¬ 
vature  (0.05  and  0.10)  collapsed  upon  one  another 
when  plotted  In  -uV/Utz  vs.  y/R  coordinates. 

They  concluded  that  cases  with  6/R  >0.05  Mre 
fundamentally  different  than  Maker  curved  cases; 
the  Mak  cases  have  evolving  shear  stress  profiles 
while  shear  stress  profiles  for  stronger  curved 
cases  have  reached  their  asymptotic  curved  state. 

The  present  study  adds  Maker  curved  heat 
transfer  and  fluid  mechanics  data,  with  recovery, 
to  the  data  base.  Since  previous  studies  attained 
various  6/R  values  by  varying  6,  the  present 
program  ms  designed  to  change  6/R  by  varying  R, 
thus  testing  the  robustness  of  the  strength  of 
curvature  parameter,  6/R. 


2.  THE  EXPERIMENTAL  APPARATUS 


The  experiment  ms  conducted  In  an  open- 
circuit,  blown-typo  wind  tunnel  constructed  with 
an  upstreaa  developing  section,  a  curved  section 
and  a  downstream  recovery  section  (Fig.  1). 

Details  of  this  wind  tunnel  are  given  In  [10]. 

The  test  channel  Is  rectangular,  68  cm  In  span  and 

11.4  cm  deep.  The  heated  test  Mil  consists  of  a 

1.4  m  long  developing  section,  a  1.4  ■  long  curved 
section  of  either  0.9  ■  or  2.1  ■  radius  of  cur¬ 
vature,  and  a  1.4  ■  long  recovery  section. 


Fig.  1.  Plan  view  of  the  curved  boundary  layer 
facility 

The  test  Mil  (convex  for  this  study)  ms 
heated  to  nominally  8°C  above  the  free-stream  air 
temperature  with  a  uniform  heat  flux  of  25S  W/m2. 
Static  pressure  on  the  test  m11  ms  uniform.  The 
three  sections  of  the  channel  wire  of  similar 
construction  except  that  the  curved  section  could 
be  bent  to  various  radii  of  curvature.  This  gave 
the  opportunity  to  systematically  Investigate  cur¬ 
vature  effects  with  a  single  Mil.  The  entry  sec¬ 
tion  provided  a  mature  turbulent  boundary  layer  In 
the  measurement  area.  Boundary  layers  entering 
the  curved  section  wre  similar  for  the  two  cases. 
Stanton  number  data  ms  corrected  for  back-side 
heat  losses,  radiation  losses,  streamwlse  conduc¬ 
tion  loss,  recovery  effects  and  the  effects  of 
temperature  and  humidity  on  fluid  properties.  The 
uncertainty  of  the  Stenton  number  data  ms  5X. 

Static  pressures  were  measured  through  0.64  mm 
diameter  taps  In  the  opposite  Mils  end  end-walls. 
In  the  curved  region,  these  pressures  Mre  used  to 
estimate  the  static  pressure  at  the  test  wall  by 
assuming  potential  flow.  Mean  velocity  profiles 
Mre  measured  with  a  0.34  mi  10  pitot  tube.  Mean 
temperature  profiles  Mre  measured' with  a  2-0 
boundary  layer  thermocouple  probe;  Turbulence 
measurements  Mre  taken  In  an  Isothermal  flow  with 
a  constant  temperature  hot-wire  anemometer;  e 
horizontal  wire  for  normal  stress  measurements  and 
a  cross-wire  for  shear  stress  measurements.  The 
digitized  anemometer  signals  Mre  computer- 
linearized  and  processed.  Signals  from  the  cross¬ 
wire  probe  Mre  digitized  slmulteneously. 

Averages  Mre  taken  over  a  period  of  40  seconds. 


2.1  Quellflcatlon  Test 

The  mean  velocity  meesured  In  the  potential 
core  of  the  developing  section,  nominally  16  m/s, 


ms  uniform  to  0.2X  and  the  mean  temperature, 
nominally  26°C  ms  uniform  to  0.03*C.  Heat  flux 
on  the  test  Mil  ms  uniform  to  within  IX  and  the 
potential  velocity  ms  uniform  to  3X  In  Cp.  Free- 
stream  turbulence  Intensity,  normalized  by  Up*  was 
0.6SX.  Spanwlse  variations  of  the  Stanton  number, 
measured  within  the  central  span  of  30  cm,  Mre 
typically  less  than  5X.  Secondary  flow  measure¬ 
ments,  taken  with  a  Conrad  probe,  shoMd  skew 
angles  of  less  than  2  degrees  within  the  curve  and 
3  degrees  within  the  recovery  section.  Energy  and 
momentum  balances  over  the  full  test  region  showed 
closure  to  within  5X  and  3%.  From  the  above,  It 
ms  felt  that  the  boundary  layer  was  sufficiently 
two-dimensional  and  that  secondary  flow  effects 
Mre  minimal. 


3.  RESULTS  AN0  DISCUSSION 

Descriptors  for  the  two  cases  discussed  herein 
are  listed  In  Table  1.  As  shown.  Cases  1  and  2 
have  similar  flow  conditions  upon  entry  to  the 
curve;  they  differ  only  In  radius  of  curvature. 


Table  1.  Descriptors  of  Cases  1  &  2  and  locations 
_ of  stat 1 ons _ 


Casi 

R 

(cm) 

too 

T.I. 

(*) 

~m - 

(0/R) 

Beginning  of 
curvature 

T7R 

(0/R) 

End  of 
curvature 

I — 

90 

18.32 

0.65 

0.03 

(0.0036) 

"OX 

(0.0054) 

2 

210 

17.12 

0.65 

0.013 

CO. 00161 

0.019 

(0.0025) 

— sTjn 

Case 

- j - 

2 

- 3 — 

- 1 - 

“5 - 

1 

*lcm) 

-29.9 

26.2 

52.8 

79.8 

Re# 

2737 

— 

3987 

4287 

4602 

2 

X(cm) 

-29.9 

0.0 

27.9 

56.1 

85.3 

R*0 

2916 

3565 

3814 

4539 

4661 

me 

Case 

s — 

— 7 - 

8 

9 

IB 

1 

XXcm) 

106.7 

133.4 

150.6 

178.6 

208.5 

*•0 

5056 

5140 

5589 

6146 

6580 

2 

X(cm) 

113.5 

141.0 

158.8 

186.7 

214.6 

Reg 

5622 

5738 

6022 

6417 

7010 

•No  data  taken  due  to  proximity  to  a  0.02S  mm 
step. 


3.1  The  Base  Case:  R  ■  90  cm;8/R  *  0.03-0.04 

The  Reynolds  number  based  on  the  momentum 
thickness  was  2740  and  the  shape  factor  was  1.42 
at  Station  1  Indicating  a  mature  turbulent  boun¬ 
dary  layer. 


Fig.  2  straws  Man  velocity  prof  list  plotted  In 
.Inner  coordinates.  Skin  friction  coefficients 
were  deduced  fro*  the  velocity  profile  (the 
Clauser  technique  [11])  and  the  Van  Driest  rela¬ 
tionship.  The  stabilizing  effect  of  convex- 
curvature  reduces  the  length  of  the  log-linear 
region  slightly  and  enhances  the  wake  greatly. 

One  sees  that  curvature  affects  only  the  outer 
region,  Y+>1S0.  The  existence  of  the  log-linear 
region  within  Y+<150  Indicates  that  the  near-wall 
layer  has  normal  structure.  This  lends  support  to 
the  use  of  the  Clauser  technique.  No  direct 
measurements  of  wall  shear  stress  were  taken.  At 
Station  10  In  the  recovery  section,  the  tog-linear 
region  appears  to  have  recovered  completely,  but 
the  wake  Is  not  growing  normally.  In  fact,  the 
wake  Is  decreasing  slightly  In  strength.  This 
Indicates  that  recovery  is  not  complete. 


Figure  2.  Mean  velocity  profiles.  Case  1,  the 
base  case. 


The  rates  of  growth  of  the  boundary  layer 
parameters,  #,  #*,  and  •,  were  reduced  by  cur¬ 
vature,  resuming  In  the  recovery  region.  The 
shape  factor  Increased  from  1.42  to  1.52  by  the 
end  of  the  curved  section,  then  decreased  slowly 
within  the  recovery  section.  The  reduced  rate  of 
growth  of  the  momentum  thickness  and  the 
Increasing  shape  factor  within  the  curve  are  asso¬ 
ciated  with  a  more  laminar-like  shape  of  the  velo¬ 
city  profiles  characterized  by  the  growth  of  the 
wake.  These  are  the  effects  of  a  decrease  In  tur¬ 
bulent  transport  In  the  outer  (wake)  regions  due 
to  curvature. 

Fig.  3  shows  mean  temperature  profiles  plotted 
In  wall  coordinates.  They  show  a  response  to  cur¬ 
vature  that  Is  similar  to  that  observed  In  the 
velocity  profiles.  The  mean  temperature  data  Is 
expected  to  follow  T*  «  Pr*Y*  near  the  wall  and 
T*  -  (Prt/0.41)1n(Y*/Ycl*)*Pr*Yei+)  In  the  log- 
linear  region.  The  mean  temperature  profiles  were 
fit  using  this  turbulent  log-linear  correlation 
assuming  Prt  and  Yr\4  are  free  parameters  to 
determine  their  values.  The  effective  average 
turbulent  Prandtl  number  deduced  from  the  profiles 
were  Increased  about  25%  by  curvature,  then 
quickly  returned  to  0.88  (Fig.  7). 


Figure  3.  Mean  temperature  profiles,  Case  1,  the 
base  case.  ' 

Fig.  4  shows  the  local  Stanton  number  and  skin 
friction  coefficient  values  comparing  the  curved- 
wall  results  to  flat-wall  values  predicted  by  the 
boundary  layer  code  STAN5  [12].  The  effects  of 
curvature  are  dramatic.  Stanton  number  responds 
Immediately  upon  entry  to  the  curve,  decreasing 
about  15%  by  x  *  30  cm  and  continuing  to  decrease 
slowly  thereafter.  The  same  trend  can  be  observed 
for  the  skin  friction  coefficient.  Recovery  of 
the  skin  friction  coefficient  Is  slower  than  that 
of  the  Stanton  number,  as  was  observed  In  [8]  and 
[9].  Stanton  numbers  Increase  to  within  13%  of 
the  flat-plate  value  80  cm  downstream  of  the  end 
of  curvature  while  skin  friction  coefficients 
remain  20%  below.  The  rapid  decrease  of  effective 
average  Prt  Is  consistent  with  this  Increase  of 
25t/Cf  In  the  recovery  section. 


Figure  4.  Stanton  number  and  skin  friction  coef¬ 
ficient  vs.  streamwlse  distance,  Case  1,  the  base 
case. 


Turbulence  Intensity  profiles  are  shown  In 
Fig.  5.  The  upstream  flat-wall  profile 
(Reg-2700)  was  found  to  lie  between  the  two  profi¬ 
les  given  by  Klebanoff  (Reg-7500),  [13,  14]  and 
compared  well  with  previous  researchers’  results. 
Between  Stations  1  and  3  the  near-wall  peak  dra¬ 
matically  falls  while  the  remainder  of  the  profile 
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(y/6  >0.1)  Is  less  affected.  The  response  to  the 
Introduction  of  convex  curvature  It  complete  by 
Station  4  and  a  self-similar  shape  continues 
through  the  reaalnder  of  the  curve.  The  o'  profi¬ 
les,  then,  Indicate  an  Immediate  response  of  tur¬ 
bulence  quantities  upon  entry  to  curvature 
followed  by  a  sore  gradual  adjustment  within  the 
curve;  this  Is  consistent  with  the  St  and  Cf/2 
trends  In  Fig.  4.  The  return  to  flat-plate  profi¬ 
les  on  the  recovery  wall  Is  different.  A  peak 
appears  at  y/6  *  0.1-0. 2  at  the  beginning  of  the 
recovery  process  and  diffuses  slowly  outward.  By 
Station  9,  the  recovery  of  turbulence  Is  complete 
beyond  y/5  ■  0.2,  but  the  near-wall  data  shows  no 
change  from  profiles  taken  at  the  beginning  of 
recovery.  The  slow  near-wall  recovery  Is  con¬ 
sistent  with  the  slow  response  of  skin  friction 
coefficient  observed  In  Fig.  4.  It  Is  expected 
that  the  profile  eventually  returns  to  the  Station 
1  shape  as  Cf/2  returns  to  flat-wall  values.  The 
recovery  length  of  the  present  experiment  was  too 
short  to  observe  this,  however. 


Figure  5.  Streamwlse-normel  turbulence  Intensity 
profiles.  Case  1,  the  base  case. 


Figure  6.'  Reynolds  shear  stress  profiles,  Case  1, 
the  base  case. 


Shear  stress  profiles  are  shown  on  Fig.  6. 

They  respond  quickly  to  curvature  but  do  not 
change  appreciably  after  Station  3. 

The  u1  profiles  also  respond  quickly  to  the  Intro¬ 
duction  of  curvature  but  continue  to  evolve  until 
Station  4.  Curvature  reduces  shear  stress  by 
about  45X  throughout  the  boundary  layer.  Recovery 
Is  slow  and  continuous  through  Station  9,  where 
the  magnitudes  of  the  shear  stress  are  comparable 
to,  or  greater  than,  the  magnitudes  at  Station  1 
for  the  outer  75X  of  the  boundary  layer  thickness. 
A  maximum  shear  stress  In  the  recovering  profiles 
occurs  away  from  the  wall,  moving  to  larger  values 
of  y/6  as  recovery  proceeds.  The  near-wall  shear 
stress  Is  affected  by  recovery  very  little.  This 
is  consistent  with  the  observed  trends  In  Cf/2. 

3.2  The  Effect  of  Radius  of  Curvature:  Intro¬ 
duction  of  a  Comparison  Case  with  R  *  210  cm; 

6/R  *  0.013-0.019 

The  comparison  case,  one  of  weaker  curvature 
than  the  base  case,  was  achieved  by  bending  the 
wall  to  a  larger  radius  of  curvature.  A  momentum 
thickness  Reynolds  number  of  2920  and  a  shape  fac¬ 
tor  of  1.42  at  Station  1  Indicate  similar  flow 
conditions  with  those  of  the  base  case.  Mean  and 
turbulence  profiles  for  Station  1  were  also  nearly 
Identical.  Stanton  number  values  on  the  pre-plate 
were  only  about  2.SX  lower  than  those  of  Case  1; 
skin  friction  values  were  about  IX  lower. 

Evolution  of  mean  velocity  and  temperature  profile 
Integral  parameters  and  shape  factor  confirm  the 
same  trends  discussed  for  the  bese  case.  A  com¬ 
parison  of  Prt  for  the  two  cates  Is  shown  on 
Fig.  7.  Less  of  en  effect  on  Pr*  for  the  weaker- 
curved  case  Is  observable.  Local  Stanton  number 
data  (Fig.  B)  show  a  1QX  reduction  from  plane-wall 
values  for  this  case  compared  to  e  20X  reduction 
for  the  base  cese  (which  has  approximately  twice 
the  value  of  6/R).  The  comparison  of  these  two 
cases  supports  the  conclusion  that  6/R  Is  an 
appropriate  scaling  parameter  over  this  range  of 
6/R. 


x(cm) 


Figure  7.  Turbulent  Prandtl  number  vs.  strea*Mlse 
distance.  Cases  1  and  2. 

The  comparison  of  turbulence  Intensity  profiles 
for  the  two  cases,  at  Station  5  (Fig.  9)  shows 
that  the  profiles  are  similar  but  that  the  Case  2 


profile  Is  still  evolving  while  the  Case  1  pro¬ 
file  hes  reached  an  asymptotic  shape  (by  Station 
4  as  shown  In  Fig.  S).  A  comparison  of  tur¬ 
bulence  Intensity  profiles  for  Station  9 


x(cm) 


Figure  8.  Stanton  number  and  shin  friction  coef¬ 
ficient  vs.  streamwlse  distance  (comparison  be¬ 
tween  Cases  1  and  2). 
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Figure  9.  Streamwlse-normal  turbulence  Intensity 
profiles  at  Station  S  (near  the  end  of  the  curved 
section),  comparison  between  Cases  1  &  2. 

(Fig.  10)  shows  that  the  Case  2  profile  recovery 
Is  more  complete;  only  the  near-wall  portion  dif¬ 
fers  from  the  Station  1  flat-wall  profiles. 
Profiles  of  shear  stress  plotted  In  -u'v'/llT’ 
vs.  y/R  coordinates  for  two  cases  are  given  on 
Fig.  11.  The  solid  line  Is  the  "asymptotic  pro¬ 
file"  (6/R  >  0.05)  suggested  by  61111s  and 
Johnston  [9],  while  the  dashed  line  Is  a  mild  cur¬ 
vature  profile  measured  by  Hoffmann  and  Bradshaw 
[15].  The  base  case  (Case  1)  results  support  the 
suggestion  by  61111s  and  Johnston  that  the 
appropriate  scaling  parameter  for  turbulent  struc¬ 
ture  for  strongly  curved  flows  Is  the  radius  of 
curvature,  rather  than  boundary  layer  thickness. 
The  base  case  data  show  that  the  flow  approaches 
this  same  asymptotic  profile  for  S/R  >  0.04. 

Though  the  base  case  profiles  appear  to  be 
approaching  an  asymptotic  shape,  It  Is  clear  that 
the  weaker-curved  Case  2  profiles  are  more  charac¬ 
teristic  of  those  observed  In  previous,  ml  Id- 


curvature  experiments.  The  two  cases,  therefore, 
represent  examples  of  weak  and  moderate-to-strong 
curvature  cases,  as  viewed  In  these  coordinates. 


y/5s*.s 


Figure  10.  Streamwlse-normal  turbulence  Intensity 
profiles  at  Station  9  (near  the  end  of  the  reco¬ 
very  section),  comparison  between  Cases  18  2- 


Flgure  11.  Turbulent  shear  stress  profiles  within 
the  curved  section,  comparison  between  Cases  1  8 
2. 


4.  SUMMARY  AND  CONCLUSIONS 

Significant  additions  to  the  data  base  are: 

1)  The  recovery  process  with  heat  transfer  was 
Investigated  for  cases  of  milder  curvature  than 
previously  studied. 

2)  The  radius  of  curvature  vis  varied  from  case 
to  case  In  the  same  study  and  data  was  taken  for 
two  6/R  values  representing  mild  and  moderately 
strong  curvature  using  the  same  apparatus. 

Conclusions  are: 

1)  Profiles  of  mean  velocity  and  temperature 
show  shortened  log-linear  regions  and  enhanced 
wake  regions  In  the  curve.  The  recovery  of  these 
mean  profiles  In  the  log-linear  region  Is  complete 
by  about  80  cm  of  recovery  length. 

2)  Curvature  Increases  turbulent  Prandtl  number 
by  about  20-25X  for  the  cases  studied. 

3)  A  rapid  decrease  of  St  and  Cf/2  at  the 
beginning  of  the  curve  Is  observed  followed  by  a 


slow  decrease  within  the  curve.  Decreases  of 
about  20  and  10*  from  axpaeted  flat-wall  valuat 
ara  obsarvad  for  Casas  1  (6/R  *  0.03)  and  2 
(6/R  »  0.013),  raspactlvaly.  Recovery  Is  extre- 
Mly  slow  with  St  recovery  leading  Cf/2  recovery, 
resulting  In  an  Increase  of  2St/Cf. 

4)  Normal  stress  profiles  show  a  fast  response 
to  the  beginning  of  curvature.  Ourlng  recovery 
an  Increase  of  turbulence  Intensity  originates 
near  the  wall  at  y/<  -  0.1-0. 2  propagating  slowly 
to  the  outer  layers. 

5)  Shear  stress  profiles  respond  quickly  to  cur¬ 
vature,  reducing  to  about  SSX  of  the  upstream 
flat-plate  values  for  S/It  *  0.03.  Recovery  occurs 
slowly,  with  the  wake  recovering  faster  than  the 
near-wall  flow. 

6)  Asymptotic  behavior  Is  observed  for  flows 
with  S/R  >0.04. 
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NOMENCLATURE 


Cf/2 


M 

R 

Re 


U 

U* 

S» 

v' 

-u'v' 

X 

Y* 

Yd* 


Skin  friction  coefficient 
Static  pressure  coefficient, 

(Psw-*’swiref)/(*rW> 

Mach  number 
Radius  of  curvature 
Reynolds  number  (using  U~,  as  charac¬ 
teristic  velocity) 

Normalized  mean  temperature, 
(VT)^CpUA 
Streamwlse  mean  velocity 
Normalized  streamwlse  velocity. 

Shear  velocity,  /r*/p 
Root -mean-square  of  f 
velocity 

Root-mean-square  of  fluctuating  velocity 

normal  to  the  test  wall 

Reynolds  shear  stress 

Streamwlse  distance 

Distance  normal  to  the  test  wall 

Normalized  distance  normal  to  the  test 

wall,  yUt/v 

Conduction  layer  thickness  In  Inner  coor¬ 
dinates 


U/UT 

fluctuating  streamwlse 


Greek  Letters 

I  Boundary  layer  thickness  based  on  99. 5%  of 

potential  flow  velocity 
6*  Displacement  thickness 

•  Momentum  thickness 


Subscripts 

p  Potential  flow 

pw  Potential  flow  value  at  the  wall 

ref  At  the  reference  station  (Station  1) 
s  Static 

sw  Static  conditions  on  the  test  wall 
w  Wall  value 

-  Free-stream  value 
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Free-Stream  Turbulence  Effects  on  Convexly  Curved 
Turbulent  Boundary  Layers 

YOU,  8.  M..  SIMON,  T.  W.,  m4  KIM.  J. 

University  of  Minnesota 
Minneapolis,  Minnesota 


ABSTRACT 

Ftee-etream  turbulence  intensity  effect*  on  the  convex-curved  tur¬ 
bulent  boundary  layer  are  investigated.  An  attached  fully  turbulent 
boundary  layer,  grown  on  a  flat  plate,  is  introduced  to  a  convex  wall  of 
constant  radius  of  curvature.  Two  cases  with  free-stream  turbulence 
intensities  of  1.85/4  and  0.6594,  taken  in  the  same  facility  with  mod¬ 
erate  curvature,  S/R  =  0.03-0.045,  are  compared.  The  two  cases  have 
similar  flow  conditions  upon  entry  to  the  curve,  thus  separating  free- 
stream  turbulence  effects  within  the  curve  from  the  other  effect*.  The 
higher  turbulence  case  shows  stronger  curvature  effects  on  C//2  and 
on  streamwise  normal  and  shear  stress  profiles  than  those  observed  in 
the  lower  turbulence  case.  The  case  with  increased  turbulence  has  a 
higher  (~  594)  C//2  value  upstream  of  the  curve  than  the  Cf/2  value 
of  the  lower  turbulence  case,  as  expected,  but  this  increase  disappears 
by  the  end  of  the  curve  where  nearly  identical  Cf/2  values  are  found. 
Similarly,  streamwise  turbulence  intensity  profiles,  differing  upstream 
of  the  curve  for  the  two  cases,  are  found  to  be  similar  near  the  end 
of  the  curve;  this  indicates  a  dominating  effect  of  curvature  over  the 
effect  of  turbulence  intensity.  Many  effects  of  curvature  observed  in 
the  lower  turbulence  intensity  case,  e.g.,  a  dramatic  response  to  the 
introduction  of  curvature,  are  also  seen  in  the  higher  turbulence  case. 

NOMENCLATURE 

Cf/2  Skin  friction  coefficient 

Cr  Static  pressure  coefficient,  2(Pti.  -  P,  i\h»U 

cr  Specific  heat 

H  Shape  factor 

P  Pressure 

Pr  Prandtl  number 

Pr,  Turbulent  Prandtl  number 

f"  Heat  flux 

R  Wall  radius  of  curvature 

Rt  Reynolds  number 

St  Stanton  number 

T  Mean  temperature 

T*  Normalised  mean  temperature,  (T.  -  7«.)lf,/(f"/pe,) 


T.L  Turbulence  intensity 

U  Streamwise  mean  velocity 

U*  Normalised  streamwies  velocity,  U/U, 

UT  Shear  velocity,  \Jrm/p 

s'  Fluctuating  component  or  root-mean  square  of  fluctuating 

streamwise  velocity 
-vfx?  Reynolds  shear  stress 

v'  Fluctuating  component  of  velocity  norms!  to  the  test  wall 

s  Streamwise  distance 

y  Distance  normal  to  the  teat  wall 

Y +  Normalised  distance  normal  to  the  test  erall, 

Conduction  layer  thickness  in  inner  coordinates 
Greek  Letters 

S  Boundary  layer  thickness  based  on  09.594  of  potential 

flow  velocity 

f  Momentum  thickness 

k  K  arm  an  constant,  0.41 

v  Kinematic  viscosity 

p  Density 

r  Shear  stress 

Subscripts 

pw  Potential  flow  value  at  the  wall 

ref  Reference 

*  Static 

w  Wall  value 

oo  Free-stream  value 

1.  INTRODUCTION 

The  present  experimental  study  was  conducted  to  enhance  the 
knowledge  of  curved  turbulent  flows  and  to  expand  the  data  base  used 
for  development  of  computational  models,  particularly  those  applied 
to  gas  turbine  design.  Despite  the  considerable  work  on  the  effects  of 
curvature,  the  present  study  is  the  first  showing  the  effect  of  turbu- 
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knot  intensity  on  convex-curved  turbulent  boundary  layer  flown.  It  is 
important  to  investigate  tbe  effect  of  turbulence  intensity  with  curva- 
ture  since  accurate  predictions  of  skin  friction  and  beat  transfer,  under 
tbe  combined  effects  of  curvature  and  free  stream  turbulence  inten¬ 
sity,  are  needed  in  tbe  design  of  fluid  machinery  where  the  blades  are 
operated  in  highly  turbulent  environments.  Pres  stream  turbulence  in¬ 
tensity  can  be  as  high  as  30%  in  tbe  gee  turbine.  The  present  study  is 
a  continuation  of  tbe  convex  curvature  turbulent  boundary  layer  study 
presented  by  You,  et  al.  [i]  in  which  two  lower  turbulence  cases  of 
different  radii  of  convex  curvature  were  tested  (S/R  —  0.013  and  0.03). 
Tbe  stronger  curved  caea  of  the  two  is  the  base  case  for  the  present 
study  which  isolates  the  turbulence  intensity  effect. 

1.1  Previous  Work 

Bradshaw  [2]  published  a  comprehensive  survey  of  the  literature 
on  the  effects  of  streamwise  curvature  in  1073  showing  that  the  ratio 
of  boundary  layer  thickness  to  radius  of  wall  curvature,  6/ R,  was  an 
appropriate  descriptor  for  the  strength  of  curvature.  Experimental 
studies  [3,4,5]  published  after  the  Bradshaw  survey  confirmed  earlier 
findings  that  streamwise  curvature  has  a  marked  effect  on  turbulence 
structure,  heat  transfer,  and  skin  friction. 

Carefully  controlled  and  detailed  experiments  on  mild  (S/R  = 
0.01)  and  strong  ( S/R  >  0.05)  convex-curved  boundary  layers  were 
performed  by  Cibeon  [0,7,8],  Simon  [9,10,11,12],  and  Gillia  [13,14]  and 
their  co-workers. 

In  1986,  You,  et  al.  [1]  added  heat  transfer  and  fluid  mechanics 
data  on  a  mildly-curved  convex  surface,  with  recovery  from  curvature, 
to  the  data  base.  This  test  was  designed  so  that  S/R  could  be  changed 
by  varying  R  in  a  single  test  facility  from  mild  (S/R  =  0.013)  to  mod¬ 
erately  strong  (S/R  =  0.03).  Though  both  cases  show  the  same  trends, 
the  effects  of  the  different  strengths  of  curvature  are  clearly  observable. 
Some  observations  were: 

(1)  The  percentage  decrease  in  St  and  Cy/2  from  expected  flat- 
wall  values  is  almost  an  order  of  magnitude  greater  than  the 
value  of  S/R]  10%  decrease  for  S/R  m  0.013  and  20%  for  S/R 
=  0.03. 

(2)  The  turbulent  Prandtl  number  deduced  from  the  mean  ve¬ 
locity  and  temperature  profiles  increases  within  the  curve  about 
20%  for  S/R  =  0.013  and  25%  for  S/R  =0.03. 

(3)  Shear  stress  profiles,  plotted  in  -u'v'/U**  vs.  y/  S  show  a 
self-similar  shape  within  the  curve,  even  for  the  weaker  curva¬ 
ture  case. 

More  details  of  recent  studies  on  the  effects  of  streamwise  curvature 
are  given  in  [15]. 

Many  investigations  have  been  reported  on  the  effect  of  free-stream 
turbulence;  some  of  the  recent  ones  are  Meier  and  Kreplin  [16],  Si- 
monich  and  Bradshaw  [17]  and  Blair  [18,19].  Meier  and  Kreplin  [16] 
concluded  that  increased  free-stream  turbulence  levels  increase  skin 
friction  and  boundary  layer  growth  rate  and  that  mean  velocity  pro¬ 
files  become  similar  in  shape  to  those  of  pipe  flow.  They  also  suggested 
that,  for  turbulence  levels  less  than  1%,  Cy/2  changes  proportionally 
to  (u'/C0L  •n<*  f°r  turbulence  levels  higher  than  1%,  Cy/2  changes 
proportionally  to  ( u'/U)m  . 

Simonich  and  Bradshaw  [17]  investigated  the  effect  of  free-stream 
turbulence  on  heat  transfer.  Their  measurements,  with  turbulence 
intensities  up  to  7%,  showed  that  Stanton  numbers  increased  with 
increasing  free-stream  turbulence. 

In  1983,  several  tests  were  conducted  by  Blair  [18,19]  in  a  tero 
pressure  gradient  two-dimensional  channel,  with  turbulence  intensities 
varying  from  0.25%  to  7%,  showing  the  effects  of  free-stream  turbu¬ 
lence  on  turbulent  boundary  layer  heat  transfer  and  hydrodynamics. 
He  concluded  that  a  free-stream  turbulence  intensity  of  6%  increased 


skin  friction  and  hast  transfer  rates  approximately  14%  and  18%,  re¬ 
spectively.  He  also  found  that  the  Reynolds  analogy  factor  (2 St/Ct) 
increased  by  just  over  1%  for  each  1%  increase  in  free-stream  turbu¬ 
lence  intensity. 

2.  EXPERIMENTAL  APPARATUS 

The  present  experiment  wss  conducted  in  an  open-circuit,  blown- 
type  wind  tunnel  constructed  with  an  upstream  developing  section,  a 
curved  section  and  a  downstream  straight  section  (Pig.  1).  Details  of 
this  wind  tunnel  are  given  in  [15].  The  test  channel  is  rectangular,  68 
cm  in  span  and  11.4  cm  deep.  The  heated  test  wall  consists  of  a  1.4 
m  long  developing  section  and  a  1.4  m  long  curved  section  of  0.9  m 
radius  of  curvature  followed  by  a  straight  section. 


Fig.  1  Plan  View  of  the  Curved  Boundary  Layer  Facility. 


The  laminar  boundary  layer  was  tripped  in  the  lower-turbulence 
case,  with  a  1.0  mm  high,  12.7  mm  wide  strip  beginning  10  cm  down¬ 
stream  of  a  suction  slot  —  a  mature  turbulent  boundary  layer  was 
established  in  the  measurement  area.  Higher  turbulence  levels  were 
obtained  by  inserting  a  coarse  grid  constructed  of  2.5  cm  aluminum 
stripe  in  a  square  array  on  10  cm  centers  at  the  entrance  of  the  nozzle. 
A  turbulence  intensity,  normalized  by  ffrw.ofl  .85%  wee  achieved  in 
the  test  region.  In  tbe  higher  turbulence  case,  the  trip  wen  removed 
and  the  boundary  layer  was  allowed  to  pass  naturally  through  transi¬ 
tion.  Stanton  numbers  were  spsnwise  uniform  to  within  4%  upstream 
of  tbe  curve  in  both  cases.  It  was  previously  shown  that  span  wise 
uniformity  is  maintained  during  natural  transition  with  2%  turbulence 
intensity  [20,21]. 

The  test  wall  was  heated  to  nominally  8®C  above  the  free  stream 
temperature  with  a  uniform  heat  flux  of  255  W/m*.  Static  pressure 
on  the  test  wall  was  uniform.  Stanton  number  data  was  corrected  for 
back-side  heat  loss,  radiation  loss,  ztrvanrwise  conduction  loaa,  recovery 
effects  and  the  effects  of  temperature  and  humidity  on  fluid  properties. 
The  uncertainty  in  the  Stanton  number  data  was  approximately  5%. 

Static  pressures  were  measured  through  0.64  mm  diameter  taps 
in  tbe  opposite  wall  and  end  walls.  In  the  curved  region,  these  pres¬ 
sures  were  used  to  estimate  tbe  static  pressure  at  the  test  wall  assum¬ 
ing  potential  flow.  Mean  velocity  profiles  were  measured  with  a  0.34 
mm  ID  pitot  tube.  Mean  temperature  profiles  were  measured  with  a 
2-D  boundary  layer  thermocouple  probe.  Turbulence  measurements 
were  taken  in  an  isothermal  flow  with  constant  temperature  hot-wire 
anemometry;  a  horizontal  wire  for  normal  stress  measurements  and 
a  cross-wire  for  shear  stress  measurements.  Tbe  digitized  anemome¬ 
ter  signals  were  computer-linearized  and  processed.  Signals  from  the 
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Fig.  4  Compariaon  of  Moan  Velocity  Profiler  at  Station  1 
(Caaaa  1  and  2). 
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Stanton  number*  and  akin  friction  coefficients  are  compared  in 
Fig.  9.  St  and  C//2  for  the  higher  turbulence  case  (T.I.  =  1.85%) 
measured  at  station  1  are  about  5%  larger  than  those  of  case  1.  Cur¬ 
vature  is  seen  to  affect  the  higher  free-stream  turbulence  data  more 
than  the  corresponding  lower  turbulence  data,  as  evidenced  by  the 
C//2  values;  although  C// 2  is  higher  for  the  higher  T.I.  case  at  the 
entrance  to  the  curve,  there  is  no  difference  in  the  values  by  the  end 
of  the  curve.  Stanton  numbers  for  both  cases  showed  similar  reduc¬ 
tions  in  spit*  of  the  enhanced  curvature  effect  on  C//2  for  the  higher 
turbulence  intensity  css*.  Turbulent  Prandtl  numbers  deduced  from 
the  mean  velocity  and  temperature  profiles  (Fig.  10)  support  the  same 
trend  as  observed  in  Fig.  9.  They  decrease  near  the  end  of  the  curved 
region  in  case  2,  while  in  case  1  they  increase  continuously  throughout 
the  curve.  This  is  consistent  with  the  increase  of  ISt/Cj  at  stations  5 
and  6  in  case  2  due  to  a  larger  effect  on  Cfj 2  than  on  St.  Turbulent 
Prandtl  numbers  are  found  by  forcing  the  mean  temperature  profile  to 
obey  the  thermal  law  of  the  wall.  The  curve  and  data  are  matched  by 
choosing  the  appropriate  values  of  Pr,  and  Y£ ,,  as  discussed  in  [15]. 
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Fig.  7  Comparison  of  Mean  Velocity  Profiles  at  Station  6 
(Cases  1  and  2). 
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Fig.  8  Shape  Factor  vs.  Streamed**  Distance  (Casas  1  and 
2)- 
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Fig.  9  Stanton  Number  and  Skin  FViction  Coefficient  vs. 
Streamwise  Distance,  Comparison  Between  Cases  1  and  2. 
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Fig.  10  1\irbul*nt  Prandtl  Number  vs.  Streamwise  Distance 
(Cases  1  and  2). 


Streamwise-normal  turbulenca  intensity  profiles  in  case  2  are  shown 
in  Fig.  It.  The  response  to  the  introduction  of  convex  curvature  shows 
the  —  trends  as  in  the  base  case.  The  reduction  of  turbulence  in 
the  in—  half  of  the  boundary  layer  from  the  flat  plate  profile  to  the 
curved  self-similar  profile  for  case  2  is  about  twice  that  of  case  1  (Fig. 
12).  Even  though  the  shapes  of  the  turbulence  intensity  profiles  at  sta¬ 
tion  1  are  different,  the  curved  self-similar  profiles  for  cases  1  and  2  are 
nearly  equal;  this  shows  that  the  stabilising  curvature  effect  dominates 
the  free-stream  turbulence  intensity  effect  (see  Fig.  12). 
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the  boundary  layer,  the  Reynolds  shear  stress  reverses  in  sign  down¬ 
stream  of  station  3 — this  was  not  observed  in  the  lower  turbulence  case 
(ease  1).  Similar  profiles,  including  reversal  of  sign,  were  previously 
reported  for  more  strongly  curved  cases  [14].  It  is  interesting  that 
there  is  a  self-similar  profile  for  y/6  <  0.5  at  stations  4  and  5  in  this 
higher  turbulence  intensity  case  (Fig.  13),  while  the  outer  part  of  the 
boundary  layer  shows  a  reversal  in  sign  of  u'v*.  In  Fig.  14,  the  shear 
stress  data  (-« 7P/U*  vs.  y/R)  approaches  the  “asymptotic*  line  for 
strong  curvature  suggested  by  Gillie  and  Johnston  [14],  as  6/R  grows. 
The  shear  stress  profile  at  6/R  -  0.042  (station  6)  is  doae  but  not 
on  this  line.  The  data  also  closely  approached  this  line  in  case  1  (low 
free-stream  turbulence  intensity)  when  6/R»»  0.04  [l]. 


Fig.  11  Streamwise-normal  Turbulence  Intensity  Profiles, 
Case  2  (T.I.=1.85%). 


Fig.  13  Reynolds  Shear  Stress  Profiles,  Case  2  (T.I.=1.85%). 


y/5 


Fig.  12  Streamwise-normal  Turbulence  Intensity  Profiles  at 
Stations  1  and  5,  Comparison  Between  Cases  1  and  2. 


Fig.  14  Turbulent  Shear  Stress  Profiles  within  the  Curved 
Section,  Case  2  (T.I.=1.85%). 


Reynolds  shear  stress  data  is  plotted  in  Fig.  13.  The  dramatic 
response  to  the  introduction  of  curvature  and  the  appearance  of  an 
asymptotic  profile  are  similar  to  those  observed  in  case  1.  The  shear 
stress  profile  at  station  1  in  case  2  shows  an  extended  “tail*  —  tur¬ 
bulent  shear  stress  extends  well  beyond  6.  The  effect  of  curvature  on 
shear  stress  profiles  within  the  curve  is  dramatic.  In  the  outer  30%  of 


4.  SUMMARY  AND  CONCLUSIONS 

Two  different  free-stream  turbulence  intensity  cases  were  studied 
to  investigate  the  effect  of  free-stream  turbulence  intensity  on  boundary 
layers  over  convex-curved  surfaces. 
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The  main  cooclnwc—  about  tba  effect  of  fraa  at  ream  turbulence 
an: 

(1)  Tba  camtan  affact  dominatae  tba  fraa  atraam  turbulence 
affact  lor  tba  caaaa  studied. 

(2)  Flat  plata  data  upstream  of  tba  earn  aritk  1.8691  tarbolaaea 
intensity  aboar  reduced  waka  atnnfth  and  about  a  5%  increase 
in  St  and  C//2  compand  to  tba  loarar  turbulence  intensity  caaa 
(0.66%). 

(3)  Tba  akin  friction  coefficient  dacraaaaa  more  for  tba  higher 
turbulanca  caaa  than  for  the  loarar  turbulence  caaa  within  the 
convex-curved  aaction.  This  relative  change  waa  not  ao  obvioua 
with  tba  Stanton  number;  thia  ia  conaiatant  with  the  loarar  Prt 
of  the  higher  turbulanca  caaa  than  of  the  lower  turbulanca  caaa 
near  the  and  of  the  curved  aaction. 

(4)  Profilaa  of  turbulanca  intenaitjr  ahow  the  tame  tranda  for 
the  taro  caaaa:  dramatic  raaponaa  to  the  introduction  of  curva¬ 
ture  and  tha  appearance  of  asymptotic  turbulent  ahaar  atraaa 
profilaa. 

(6)  Similar  atreamariae  normal  turbulanca  intenaitjr  profilaa  for 
tha  taro  caaaa  ware  found  naar  tba  and  of  tha  curve  in  apite  of 
the  different  profilaa  u  pat  ream  of  tha  curve. 

(6)  Rejrnolda  ahaar  atraaa  profilaa  for  tha  higher  turbulanca  caaa 
an  dramatically  influenced  by  curvature  —  ahaar  atraaa  ravaraaa 
aign  for  y/6  >  0.7  within  tha  curve.  Thia  waa  alao  reported  ia 
etrongly-curved  convax-waU  caaaa  by  Gillia  and  Johnaton  [14]. 
A  aelf-eimilar  ahaar  atraaa  profile  ia  -iiV/U*^  va.  y/6  coordi- 
natac  on  the  curved  wall  waa  realised  for  both  caaaa  inaide  50% 
of  tha  boundary  layer  thichneae. 
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ABSTRACT 

The  effects  of  streamwise  convex  curvature,  recovery,  and  free- 
itream  turbulence  intensity  on  the  turbulent  transport  of  heat 
and  momentum  in  a  mature  turbulent  boundary  layer  is  inves¬ 
tigated.  A  special  three-wire  hot-wire  probe  developed  for  ti  lls 
purpose  k  described.  Two  cases  with  free-stream  turbulence 
levels  of  O.M%  and  2.0%,  taken  in  the  same  facility  with  moder¬ 
ate  strength  of  curvature,  S/R  m  0.03,  are  compared.  Profiles 
of  u7P,  t',  TP,  and  Pi7  are  dramatically  reduced  within  the 
curve,  with  asymptotic  profiles  being  achieved  quickly  for  the 
low  T.l.  case.  Recovery  occurs  rapidly,  with  the  profiles  often 
overshooting  fiat-wall  upstream  values.  Increased  free-stream 
turbulence  has  the  effect  of  increasing  the  profiles  through¬ 
out  the  boundary  layer  on  the  flat  developing  wall.  Profiles 
agreeing  with  the  asymptotic  profiles  of  the  low  T.I.  case  are 
observed  by  the  end  of  the  curve,  however,  illustrating  the 
dominance  of  curvature  over  free-stream  turbulence  intensity. 
For  the  higher  T.I.  case,  a  reversal  in  the  sign  of  vV  in  the 
outer  half  of  the  boundary  layer  k  observed,  leading  to  neg¬ 
ative  values  of  the  turbulent  Prandtl  number  in  this  region. 
Thk  indicates  a  breakdown  in  Reynolds  analogy. 

NOMENCLATURE 

I  offset 

4  diameter  of  wire 

H  shape  factor 

k  conductivity  of  air  or  correction  coefficient  of 

Champagne 

I  active  length  of  wire 

us  elope 

No  Nuseelt  number 

OB  overheat  ratio  of  wire  based  on  free-stream 

temperature 


P  Static  pressure 
Prt  turbulent  Prandtl  number 

R  radius  of  curvature 

Re  Reynolds  number 

St  Stanton  number 

(  instantaneous  temperature 

T  time-averaged  temperature 

TI  turbulence  intensity 

u  instantaneous  streamwise  velocity 

V  time-averaged  streamwise  velocity 

v  instantaneous  cross-stream  velocity 

V  Voltage 

v>  instantaneous  cross-span  velocity 

x  streamwise  distance 

y  distance  normal  to  wall 

a  angle  between  main  flow  direction  and  direction 

normal  to  wires 

6  boundary  layer  thickness  based  on  09.5%  of 
free-stream  velocity 
p  density 

p  dynamic  viscosity 

Subscripts 

sne  anemometer 

e//  effective 

pw  potential  value  at  wall — reference  value 

tw  static  value  at  the  wall 

oo  free-stream  value 

$  momentum  thickness 

u>  wire  or  wall  value 

Superscripts 

i  fluctuating  component  or  rat,  depending  on 

context 

«  per  unit  area 


*  per  unit  time  (overdot) 

—  time  average  (over  line) 

INTRODUCTION 

Measurement  Technique! 

The  follow  inf  describe*  the  development  of  a  three-element 
hot-wire  probe  to  measure  the  fluctuating  component*  of  ve¬ 
locity  and  temperature  in  2-D  boundary  layers,  and  it*  use  in 
a  turbulent  boundary  layer  over  a  convex  surface.  The  probe 
is  based  on  a  design  by  Blair  and  Bennett  (1984).  These  mea¬ 
surements  are  important  to  the  gas  turbine  industry  for  they 
support  the  development  of  turbulent  transport  models  used 
to  predict  the  convective  heat  transfer  to  the  turbine  blades. 
These  measurements  are  difficult  and  the  data  base  is  small. 
In  fact,  the  effects  of  moderate-to-strong  curvature,  recovery, 
and  turbulence  intensity  on  the  turbulent  transport  of  haat  has 
been  studied  here  for  the  first  time. 

Previous  investigators  have  measured  fluctuating  velocity 
and  temperature.  Representative  examples  of  their  work  are 
reviewed  below:  A  description  of  the  use  of  multiple  overheats 
in  hot-wire  anemometry  to  separate  the  temperature  and  ve¬ 
locity  components  in  a  flow  was  given  by  Corrsin  (1947).  Be 
measured  temperature  fluctuations  in  the  mixing  of  heated  gas 
streams  and  concentration  fluctuations  in  the  mixing  of  differ¬ 
ent  gases  at  constant  temperature.  His  use  of  a  single  wire, 
however,  precluded  the  possibility  of  simultaneous  tempera¬ 
ture/concentration  and  temperature/velocity  measurements. 
Sakao  (1973)  used  two  parallel  Sum  wires  0.2  mm  apart  op¬ 
erating  at  constant  but  different  temperatures  to  simultane¬ 
ously  measure  instantaneous  streamwise  velocity  and  temper¬ 
ature.  A  good  response  to  frequencies  to  300  H*  was  reported. 
Hishida  and  Nagano  (1978)  used  the  same  configuration  but 
operated  the  front  wire  at  constant  current  and  the  rear  wire 
at  constant  temperature.  No  contamination  from  the  front 
wire  was  reported  and  the  frequency  response  was  reported  to 
be  6  kHs. 

Chen  and  Blackwelder  (1978)  used  a  triple  wire  probe  to 
measure  two  components  of  velocity,  and  temperature.  The 
probe  consisted  of  a  convent''  *1  cross-wire  with  a  resistance 
thermometer  placed  directly  in  front  of  the  center  of  the  cross¬ 
wires  and  in  the  plane  perpendicular  to  them.  The  probe  di¬ 
mensions  and  frequency  response  were  not  presented.  Smils 
and  Perry  (1981)  used  a  technique  similar  to  that  used  by 
Chen  and  Blackwelder  (1978),  but  with  the  temperature  sens¬ 
ing  element  placed  in  a  plane  next  to  and  parallel  with  the 
cross  wires.  The  methodology  for  separating  the  velocity  and 
temperature  components  was  presented,  but  no  measurements 
were  taken.  A  third  configuration,  used  by  Gibson  and  Ver- 
riopoulos  (1984),  consisted  of  a  conventional  cross-wire  with  a 
third  "cold”  wire  positioned  between  and  in  a  plane  parallel  to 
the  planes  of  the  cross-wires.  A  good  response  to  frequencies 
to  6.8  kHz  was  reported. 

A  4-wire  probe  designed  to  measure  three  components 
of  velocity  and  temperature  was  described  by  Fabris  (1978). 
Three  wires  were  operated  in  the  constant  temperature  mode, 
while  the  fourth  was  of  constant  current.  All  wires  were  0.625 


nm  in  diameter,  minimizing  contamination  between  the  wires. 
The  velocity  components  and  temperature  were  found  by  si¬ 
multaneously  solving  the  four  non-linear  response  equations 
for  the  sensors  on  a  computer.  A  good  response  to  frequencies 
to  2000  Hz  was  reported. 

Blair  and  Bennett  (1984)  described  a  3-wire  probe  for  use 
in  2-D  boundary  layers.  The  probe  consists  of  three  wires  lo¬ 
cated  in  three  parallel  planes.  The  two  outer  wires  are  orthog¬ 
onal  forming  an  X-array,  while  the  third  center  wire  is  parallel 
to  one  of  them.  AH  three  are  operated  at  constant  tempera¬ 
ture  with  the  center  wire  operated  at  a  much  lower  tempera¬ 
ture  than  the  two  outer  wires.  The  instantaneous  temperature 
is  found  from  the  two  parallel  wires,  and  the  instantaneous 
velocity  from  the  two  orthogonal  wires.  A  disadvantage  to 
this  design  is  that  the  measurement  of  temperature  is  not  di¬ 
rect.  The  temperature  must  be  determined  by  simultaneously 
solving  two  closely  coupled  non-linear  equations.  The  advan¬ 
tages  are  numerous,  however,  for  only  constant  temperature 
anemometers  are  used  and  frail  sub-micron  diameter  constant 
current  wire*  are  avoided.  A  good  response  to  frequencies  to 
80  kHz  was  reported.  For  this  reason,  the  scheme  of  Blair  and 
Bennett  (1984)  was  chosen  for  the  present  study.  The  wire  con¬ 
figuration  of  the  probe  used  in  this  study  is  shown  on  Fig.  1. 


Mtwurtatnn 

The  probe  was  used  to  measure  profiles  of  turbulence  quanti¬ 
ties  in  a  sero  pressure  gradient,  turbulent  boundary  layer  in¬ 
fluenced  by  convex  curvature.  Streamwise  curvature  is  known 
to  markedly  affect  the  structure  of  turbulent  boundary  layers. 


Fig.  1  ~  Schematic  of  Wire  Configuration 


The  changes  in  boundary  layer  turbulence  intensity,  wall  skin 
friction,  and  Stanton  number  from  the  corresponding  fiat-wall 
values  are  generally  an  order  of  magnitude  greater  than  the 
magnitude  of  the  strength  of  curvature,  6/R.  A  number  of 
studies  concerning  the  effects  of  curvature  on  the  mean  flow 
and  surface  heat  flux  are  described  by  You,  Simon  and  Kim 
(1986a).  Early  work  consists  almost  entirely  of  mean  temper¬ 
ature  profiles  and  surface  heat  flux  measurements.  Simon  and 
Moffatt  (1982a,  1982b)  were  the  first  to  report  heat  transfer 
data  in  which  the  details  of  the  velocity  field  were  know  n- 
given  by  Giliis  and  Johnston  (1980).  They  were  also  the  first 
to  document  the  recovery  from  curvature.  Their  strength  of 
curvature  was  l/R  *  0.1,  indicating  strong  curvature-  You,  Si¬ 
mon  and  Kim  (1986a)  added  heat  transfer  and  fluid  mechanics 
data  on  a  mildly-curved  convex  surface,  with  recovery,  to  the 
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database.  Two  strengths  of  curvature  (6/R  =  0.013,0.03),  ob¬ 
tained  by  bending  the  test  wall,  were  investigated.  This  study 
by  You,  Simon,  and  Kim  (1986a)  is  the  basis  for  the  present 
study.  Some  of  their  results  were:  1).  Profiles  of  the  mean  ve¬ 
locity  and  temperature  showed  shortened  log-linear  regions  and 
enhanced  wake  regions  in  the  curve.  2).  Curvature  increased 
the  turbulent  Prandtl  number,  deduced  from  the  mean  veloc¬ 
ity  and  temperature  profiles,  by  about  20-25%.  3).  A  rapid 
decrease  in  St  and  C//2  at  the  beginning  of  the  curve  was 
observed  followed  by  a  slow  decrease  within  the  curve.  De¬ 
creases  of  about  20  and  10%  from  expected  fiat-wall  values 
were  observed  for  5/R= 0.03  and  0.13,  respectively.  Recovery 
was  extremely  slow,  with  Si  recovery  leading  C//2  recovery. 
4).  Normal  stress  profiles  responded  rapidly  to  the  beginning  of 
curvature.  During  recovery,  an  increase  in  turbulence  intensity 
originated  near  the  wall  at  y/6  m  0.15  propagating  slowly  to 
the  outer  layers.  5).  Shear  stress  profiles  responded  quickly  to 
curvature,  reducing  to  about  55%  of  the  upstream  flat-wall  val¬ 
ues  for  6/R  =  0.03.  Recovery  occurred  slowly,  with  the  wake 
recovering  faster  than  the  near-wall  flow.  Asymptotic  shear 
stress  behavior  (first  observed  by  GUlis  and  Johnston — 1980) 
was  found  when  6/R  >  0.04. 

There  are  very  few  studies  in  which  the  effect  of  free- 
stream  turbulence  intensity  on  boundary  layer  heat  transfer 
has  been  investigated.  The  study  by  You,  Simon  and  Kim 
(1986b)  was  the  first  and  only  work  to  document  free-stream 
turbulence  effects  (TI=0.68%  and  1.85%)  on  heated  curved 
flows.  No  recovery  data  was'  taken  due  to  the  disappearance 
of  the  potential  core  in  recovery  for  the  higher  turbulence 
case.  This  turbulence  intensity  study  by  You,  Simon  and  Kim 
(1986b)  also  provides  a  base  for  the  present  measurements. 
Some  results  were:  1).  The  curvature  effect  dominated  the 
turbulence  intensity  effect.  2).  The  skin  friction  coefficient,  de¬ 
duced  by  use  of  the  Clauser  technique,  decreased  more  within 
the  curve  than  did  Stanton  number.  3).  Similar  streamwise 
normal  turbulence  intensity  profiles  for  the  two  cases  of  dif¬ 
ferent  TI  were  found  near  the  end  of  the  curve  in  spite  of  the 
different  profiles  upstream  of  the  curve.  4).  Profiles  of  shear 
stress  showed  the  same  trend  for  both  cases:  a  dramatic  re¬ 
sponse  to  the  introduction  of  curvature  and  asymptotic  turbu¬ 
lent  shear  stress  profiles  by  the  end  of  the  curve.  A  reversal  in 
the  sign  of  the  shear  stress  was  seen  in  the  high  TI  case  at  the 
entry  to  the  curve. 

Gibson  and  Verriopoulos  (1984)  were  the  first  and,  until 
the  present  study,  the  only  workers  to  take  turbulence  measure¬ 
ments  in  a  heated  curved  (6/R  =  0.01)  flow.  Their  measure¬ 
ments  indicated  that  v't'  was  affected  by  curvature  more  than 
was  uV.  An  initial  sharp  fall,  then  a  continuous  decrease  of 
v't '  was  observed.  I7!7  was  less  strongly  affected  by  curvature. 
It  responded  more  slowly  to  the  step  change  in  curvature.  Its 
profiles  at  the  last  two  stations  in  the  curve  were  very  nearly 
similar.  Pr,  values  were  scattered,  but  the  authors  felt  that  a 
rise  in  Prt  was  observable.  This  is  consistent  with  the  values 
deduced  from  mean  velocity  and  temperature  profiles,  and  with 
the  results  of  Simon  and  Moffatt  (1982a,  1982b),  and  You,  Si¬ 
mon,  and  Kim  (1986a).  A  final  conclusion  was  that  turbulence 


in  the  boundary  layer  is  modified  by  wall  curvature  such  that 
heat  transfer  was  stabilised  more  effectively  than  momentum 
transfer. 


PROBE  DESIGN  AND  DEVELOPMENT 

In  designing  special  purpose  hot-wire  arrays,  consideration  must 
be  given  to  a  multitude  of  factors,  some  of  which  are  spatial 
averaging,  prong  and  shaft  interference,  end  conduction,  sen¬ 
sor  cross-talk,  and  survivability.  The  probe  for  the  present 
study  was  constructed  using  the  guidelines  presented  by  Blair 
and  Bennett  (1984).  The  wires  are  2.5  nm  diameter  platinum- 
plated  tungsten  with  an  active  length-to-diameter  ratio  of  200 
and  a  separation  distance  of  0.35  mm  (Fig.  1).  The  ends  of  the 
wires  are  plated  to  reduce  prong  interference  and  end  conduc¬ 
tion  loss. 

The  response  equation  of  each  sensor  is  assumed  to  be  of 
the  form 


Nu  *  Ax  +  BiRt0  43t 

which  is  a  slight  variation  on  King’s  law.  Substituting  in  the 
definition  of  JVu  and  Rt,  modelling  the  property  variations  as 

and  re-arranging  leads  to  the  sensor  response  equation 


=  At™6  + 


&  yt 
(tu,  -  «oo)  *"* 


where  A  and  B  are  assumed  constant.  This  is  the  equation 
used  by  Blair  and  Bennett  (1984)  in  modelling  sensor  response. 
The  authors  have  found,  however,  that  A  and  B  are  slight 
functions  of  temperature.  This  variation  is  incorporated  into 
the  response  equation  using  a  least-squares  fit  to  the  calibration 
data.  The  final  response  equation  is  given  by 


u0.«ss 


(rrtjit  oo  +  4a)<2o76  + 


(mat  OP  ^B)  yj 
(K  -  too)  *n‘ 


The  two  assumptions  made  when  reducing  the  data  are  1).  The 
boundary  layer  is  two-dimensional,  and  2).  The  instantaneous 
velocities  seen  by  the  two  parallel  wires  are  equal.  The  first 
assumption  is  needed  since  the  third  component  of  velocity  (v) 
cannot  be  determined.  If  the  probe  is  aligned  with  the  flow, 
however,  and  if  the  boundary  layer  is  two  dimensional,  the 
tvcomponent  makes  only  second-order  contributions  to  sensor 
response  and  may  be  safely  neglected.  If  the  second  assumption 
holds,  the  velocity  term  («,//)  may  be  eliminated  from  the 
response  equations  for  the  two  parallel  wires,  resulting  in  an 
equation  in  which  the  sole  variable  is  the  ambient  temperature 
too-  The  ambient  temperature  is  found  by  iteratively  solving 
for  too  using  the  Newton-Raphson  method.  The  instantaneous 
velocities  a  and  v  may  then  be  found  from  the  signals  of  the 
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two  oat «r  orthogonal  wire*  using  Champagne’s  form  of  the  k * 
relations: 

j  m  («  eo*a  +  v  etna)*  +  k*(u  sina  -  v  eo*a)* 


a*^  «  («  eosa  -  v  lino)*  +  k*(u  sina  -  v  eosa)* 

Knowing  the  instantaneous  values  of  u,  v,  and  t,  the  rms  fluc¬ 
tuation  quantities  («',  v',  and  t')  and  their  cross  correlations 
(*V,  ult\  and  Tt1)  may  be  determined.  The  probe  was  cali¬ 
brated  as  a  function  of  both  velocity  and  temperature. 

Qualification  of  the  probe  was  in  a  sero-  pressure  gradient 
flat-plate  turbulent  boundary  layer  with  a  momentum  thick¬ 
ness  Reynolds  number  of  3670  and  a  uniform  wall  heat  flux 
of  160  W/m*.  The  boundary  layer  thickness  and  free-stream 
velocity  were  3.35  cm  and  15.5  m/s,  respectively.  The  probe 
was  traversed  across  the  boundary  layer  and  measurements  of 
*V,  I',  I1?,  and  v't'  were  made.  These  quantities,  normalised 
on  the  free-stream  velocity  ((/«,)  and  the  wall  to  free-stream 
temperature  difference  {Tw  -  Too),  were  compared  with  the 
boundary  layer  data  of  Blair  and  Bennett  (1984)  and  Gibson 
and  Verriopouloe  (1964).  The  three  data  sets  agreed  well  ex¬ 
cept  in  the  vicinity  of  the  wall  (y/6  <  0.25)  where  the  present 
data  increased  beyond  those  of  other  researchers.  Insight  into 
this  discrepancy  may  be  found  from  measurements  of  appar¬ 
ent  t'.  Apparent  t'  is  the  t'  measured  by  the  probe  in  an  un¬ 
it  ea  ted  boundary  layer,  and  is  a  result  of  the  two  parallel  wires 
not  experiencing  the  same  velocity.  If  the  probe  were  perfect, 
the  apparent  V  across  the  boundary  layer  would,  of  course,  be 
sero.  A  large  apparent  t'  would  indicate  serious  problems  with 
the  probe.  Results  of  the  measurements  for  Tl=0.68%  and 
TI«2.0%,  shown  on  Fig.  3,  show  that  the  apparent  t'  values 
are  a  small  percentage  of  the  actual  values  in  the  outer  part  of 
the  boundary  layer,  but  rise  rapidly  as  the  wall  is  approached. 
This  is  expected  since  the  eddy  size  decreases  with  decreasing 
normal  distance  from  the  wall,  increasing  the  likelihood  of  the 
two  parallel  wires  not  seeing  the  same  velocity.  It  may  also  be 
seen  that  turbulence  intensity  has  little  effect  on  apparent  t’, 
further  suggesting  that  distance  from  the  wall  is  the  controlling 
parameter.  The  effect  of  apparent  t'  on  the  actual  t'  is  smaller 
than  the  curves  suggest  since  the  apparent  t'  affects  the  actual 
data  in  a  root-sum-square  manner. 

The  turbulent  Prandtl  number  (Pr«)  deduced  from  the 
measurements  of  uV  and  v’t7  were  compared  with  the  data  of 
Blair  and  Bennett  (1984)  and  Gibson  and  Verriopouloe  (1984). 
There  was  some  unavoidable  scatter,  but  the  three  data  sets 
agreed  within  their  uncertainties  (*  20%). 

The  effect  of  overheat  ratio  (OH)  on  apparent  t',  presented 
in  Fig.  3,  shows  that  increasing  OH  increases  Furthermore, 
very  low  OH’*  are  needed  to  achieve  low  apparent  t1.  Why  this 
is  the  case  is  not  presently  known,  but  it  is  thought  to  be  tied  to 
the  disappearance  of  the  velocity  dependence  and  the  reduction 
of  the  solution  matrix  stiffness  as  the  wire  temperature  1*  re¬ 
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Fig  2  -  Comparison  o)  apparent  r  with  the  1  m  a  heated  bojnoa-,  laye- 


Fig  3  -  Eflect  ot  oeertieat  rat®  (OH)  on  apparent  t' 


duced.  There  was  some  worry  that  the  low-  frequency  response 
of  the  probe  would  affect  the  measurements,  but,  judging  from 
the  qualification  data,  the  probe  seems  to  follow  fluctuations 
associated  with  important  scales  of  the  flow.  A  good  response 
to  frequencies  as  large  as  a  few  hundred  hertz  is  expected. 

As  a  last  check  on  the  probe,  the  measurements  of  shear 
stress  were  taken  with  the  two  outer  orthogonal  wires  and  again 
with  the  two  inner  orthogonal  wires  to  check  for  prong  inter¬ 
ference  and  cross-talk.  The  difference  in  the  two  data  sets  was 
well  within  the  uncertainty,  and  the  curves  agreed  well  with 
the  data  of  Gibson  and  Verriopouloe  (19S4).  No  evidence  of 
interference  or  other  forms  of  contamination  was  evident. 

In  summary,  the  probe  performs  well  except  in  the  near 
vicinity  of  the  wall.  The  uncertainty  in  the  measurements  is 
estimated  at  15%  for  the  correlations  and  20%  for  Prt. 

DESCRIPTION  OF  THE  TEST  FACILITY 

The  present  experiment  was  conducted  in  an  open-circuit,  blow  n- 
type  wind  tunnel  constructed  with  an  upstream  developing  sec¬ 
tion,  a  curved  section  and  a  downstream  recovery  section.  A 


Fig  4  -  Plan  «•*  ol  Pit  curved  boundary  lay  O'  lacilily 


wu  adjusted  to  within  3%  of  the  mean.  The  reference  pressure 
was  taken  to  be  the  static  pressure  at  station  1 . 

Low  Tl  Cat*  (TU0.$B%) 
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schematic  of  the  test  facility  is  shown  on  Fig.  4.  Details  of  the 
tunnel  are  given  in  You,  Simon,  and  Kim  (1886a  and  1886b). 
The  test  channel  is  rectangular,  68  cm  in  span  and  11.4  cm 
deep.  The  heated  test  wall  consists  of  a  1.4  m  long  develop¬ 
ing  section  and  a  1.4  m  long  curved  section  of  0.8  m  radius  of 
curvature  followed  by  a  1.4  m  long  straight  recovery  section. 

The  laminar  boundary  layer  was  tripped  in  the  lower- 
turbulence  case,  with  a  1.0  mm  high,  12.7  mm  wide  strip  be¬ 
ginning  10  cm  downstream  of  a  suction  slot — a  mature  bound¬ 
ary  layer  was  established  in  the  measurement  area.  The  free¬ 
st  ream  turbulence  intensity,  normalized  on  Vfv,  was  0.68%. 
Higher  turbulence  levels  were  obtained  by  inserting  a  coarse 
grid  constructed  of  2.5  cm  aluminum  strips  in  a  square  array 
on  10  cm  centers  at  the  entrance  of  the  nozzle.  The  boundary 
layer  trip  was  removed  and  the  boundary  layer  was  allowed  to 
pass  naturally  through  transition.  A  turbulence  intensity  of 
2.0%  was  achieved  in  the  test  section  with  the  grid  in  place. 
Stanton  numbers  were  spsnwise  uniform  to  within  4%  of  the 
centerline  value  upstream  of  the  curve  in  both  cases. 

The  test  wall  was  heated  to  nominally  6CC  above  the  free- 
stream  temperature  with  a  uniform  heat  flux  of  160  W/m3. 
Static  pressure  on  the  test  wall  was  uniform.  Static  pressures 
were  measured  through  0.64  mm  diameter  taps  in  the  opposite 
wall  and  end  walls.  In  the  curved  region,  these  pressures  were 
used  to  estimate  the  static  pressure  at  the  test  wall  assuming 
potential  flow. 

RESULTS  AND  DISCUSSION 

Measurements  of  t\  v7!7,  and  xTt7  were  taken  at  the  flat 
upstream  developing  station  (station  1),  three  stations  in  the 
curve  (stations  3,  4,  and  6),  and  two  stations  in  the  recovery 
(stations  8  and  10).  Data  was  not  taken  beyond  station  6  in  the 
high  TI  case  due  to  merging  of  the  test-wall  and  opposite-wall 
boundary  layers.  Descriptors  for  the  two  cases  are  presented 
on  Table  1.  The  flexible  outer  walls  were  adjusted  such  that 
the  pressure  coefficient  (Cr)  defined  as 


High  Tl  Case  (Tl=2.0%) 
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x(cm) 
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H 

1 

-29.87 

2640 

1.35 
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26.16 

4053 

1.41 

3 

52.83 

4155 

1.42 

Table  1  -  Sum.mry  ol  boundary  layer  parameters 


Low  TI  Case  (TI  ■  0.68%) 

Shear  stress  profiles  are  shown  on  Fig.  5.  The  profiles  are  seen 
to  respond  quickly  to  curvature.  A  dramatic  reduction  at  the 
onset  of  curvature  and  the  appearance  of  what  appears  to  be 
an  asymptotic  state  in  the  curve  occurs.  Curvature  is  seen 
to  reduce  the  shear  stress  by  about  45%  below  the  fiat  wall 
values  throughout  the  curve.  Recovery  of  shear  stress  is  slow  , 
with  recovery  aeemingly  complete  in  the  wake  by  station  10 


(See  Fig.  4  for  station  locations).  You,  et  al.  (1886s)  h*ve 
shown  that  the  near-wall  values  of  u'v’  recover  very  slowly.  All 


the  trends  ere  consistent  with  the  observations  of  Gillis  and 
Johnston  (1983)  and  You,  et  al.  (1986a). 

Profiles  of  t'  are  plotted  on  Fig.  6.  The  flat  upstream 
values  (station  1)  agree  with  the  measurements  of  Blair  and 
Bennett  (1984)  and  Gibson  and  Verriopoulos  (1984)  within 
the  uncertainty  of  the  measurement,  except  near  the  wall.  The 
profiles  within  the  curve  assume  an  asymptotic  shape  for  y/t  < 
0.8.  A  slight  evolution  of  the  profiles  for  y/t  >  0.5  is  evident. 
Recovery  has  a  dramatic  effect  on  t'.  Values  at  station  10  are 
seen  to  overshoot  the  station  1  values.  It  is  believed  that  the 
profile  eventually  returns  to  the  station  1  shape.  The  recovery 
length  was  too  short  to  observe  this  in  the  present  experiment, 
however. 
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Fig  6  ~  Eltect  of  curvature  and  recovery  on  r.  Th>0  66% 


The  effects  of  curvature  and  recovery  on  the  streamwise 
turbulent  heat  flux  n't7  are  shown  on  Fig.  7.  Station  1  values 
agree  well  with  the  data  of  Gibson  and  Verriopoulos  (1984). 
The  effect  of  curvature  on  v't'  is  dramatic.  The  profiles  in  the 
curve  “snap*  into  an  asymptotic  shape  for  y/t  <  0.5  while 
slow  evolution  of  the  profiles  for  y/t  >  0.5  may  be  observed. 
Recovery  is  seen  as  an  overshoot  of  u't'  beyond  the  station  1 
values,  although  the  profile  is  expected  to  eventually  return  to 
the  station  1  values. 

Profiles  of  cross-stream  turbulent  heat  flux  v't'  are  shown 
on  Fig.  8.  The  effects  of  curvature  and  recovery  on  v't7  are  very 
similar  to  those  observed  for  v't' — a  dramatic  reduction  in  the 
curve  to  an  asymptotic  shape,  followed  by  an  overshoot  in  the 
recovery.  Values  of  an  effective,  extrapolated  v't7  at  the  wall 
as  calculated  from  the  wall  heat  flux  measurements  are  located 
on  the  y  »  0  axis  in  Fig.  8.  The  wall  values  and  the  profile 
measurements  agree  very  well  at  station  1  and  in  the  curve. 
The  values  of  v't7  in  the  recovery  region,  however,  rise  above 
the  wall  heat  flux  values,  especially  for  station  10.  This  may 
be  seen  more  clearly  on  Fig.  9  where  vT'  has  been  normalised 
on  the  wall  heat  flux.  The  profiles  at  station  1  and  within 
the  curve  approach  unity  near  the  wall  while  the  profiles  in  the 
recovery  section  rise  above  unity.  It  is  believed  that  the  profiles 
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turn  down  and  approach  unity  as  the  wall  is  approached.  The 
measurements  cannot  be  taken  this  close  to  the  wall,  however. 
A  positive  slope  in  v't'  suggests  mixing  of  increasing  intensity 
with  increasing  distance  away  from  the  wall. 


Profiles  of  turbulent  Prandtl  number  (Pr,)  ere  shown  on 
Fif.  10.  The  data  show  some  unavoidable  scstter  which  makes 
it  difficult  to  draw  firm  conclusions  about  the  effect  of  curva¬ 
ture.  Gibson  and  Verriopoulos  (1984),  You,  et  al.  (1986a),  and 
Simon  and  Moffatt  (1982a,  and  1982b)  stated  that  Prt  in  the 
log- linear  ration  rises  due  to  convex  curvature.  The  You,  et 
al.  (1986)  and  Simon  and  Moffatt  (1982a,  and  1982b)  studies 
cite  indirect  Prt  values  deduced  from  mean  temperature  and 
velocity  profiles  which  represent  aver  ate  values  for  the  turbu¬ 
lent  core.  Because  of  scatter  in  Fig.  10  the  only  conclusion 
which  can  be  drawn  from  this  data  is  that  the  pre-plate  and 
curved-flow  Prt  values  are  not  vastly  different  (>  2096)  from 
one  another  whereas  the  data  at  the  end  of  the  recovery  section 
is  lower. 
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A  plot  of  the  triple  product  uV5  representing  the  cross¬ 
stream  flux  of  turbulent  stress  is  presented  on  Fig.  11.  The 
fiat-wall  data  is  compared  to  the  data  of  Gibson  and  Verriopou¬ 
los  (1984),  and  is  seen  to  be  higher  throughout  the  boundary 
layer.  Why  this  is  the  case  is  not  presently  known.  The  effect 
of  curvature  on  uV*  is  dramatic.  Values  of  u'v'1  are  reduced 
to  1596  to  2596  of  the  flat  wall  values  by  station  6,  indicat¬ 
ing  that  the  transport  of  str-ss  from  the  near-wall  production 
layer  to  the  outer  flow  is  virtually  shut  off.  A  sharp  drop 
at  the  introduction  of  curvature  followed  by  a  slow  continued 
decrease  is  evident.  The  profiles  are  seen  to  recover  slowly. 
The  near-wall  values  remain  low  (which  is  consistent  with  the 
downturn  in  the  shear  stress  profiles  near  the  wall)  within  the 
recovery  section,  as  discussed  earlier.  The  peak  in  the  pro¬ 
files  increases  and  is  seen  to  move  away  from  the  wall  with 
increasing  streamwiae  distance.  The  profile  appears  to  become 
negative  for  y/6  <  0.35,  indicating  diffusion  of  stress  toward 
the  wall. 

Profiles  of  the  triple  product  v^t\  the  cross-stream  flux 
of  Vt1,  is  presented  on  Fig.  12.  The  flat-wall  values  are  seen  to 
be  somewhat  higher  than  the  data  of  Gibson  and  Verriopoulos 
(1984),  as  was  the  case  for  «V*.  Curvature  is  seen  to  affect 
v'V  in  much  the  same  way  as  with  u'v'*.  A  sharp  reduction 


Fig  12  E«eci  of  curvature  and  recovery  on  v‘i  Ti.O  66' . 


in  v'2t'  followed  by  a  slow  evolution  to  values  15  to  2o9c  the 
magnitude  of  the  upstream  flat-wall  values  is  evident.  Recovery 
of  the  profiles  with  the  peak  moving  away  from  the  wall  is  again 
seen,  and  a  reversal  in  sign  of  v'2t'  is  seen  for  y/6  <  0.35,  as 
was  seen  for  uVJ.  The  latter  indicates  that  the  turbulent  heat 
flux  decreases  as  the  wall  is  approached.  This  is  consistent 
with  the  observed  downturn  in  tPP  upon  recovery. 

higher  free-stream  turbulence  case 
are  presented  on  Fig.  13.  Higher  turbulence  intensity  has  the 
effect  of  increasing  the  shear  stress  throughout  the  boundary 
layer.  This  is  expected.  The  profiles  immediately  assume  an 
asymptotic  shape  within  the  curve  for  y/6  <  0.5,  as  in  the 
low  TI  case.  The  asymptotic  profile  for  the  low  T1  case  is  also 
plotted  on  Fig.  13  where  it  is  seen  to  agree  well  with  the  high 
TI  curved  flow  data.  A  reversal  in  the  sign  of  the  shear  stress 
is  seen  to  occur  in  the  outer  part  of  the  boundary  layer.  This 
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reversal  in  sign  has  also  been  observed  by  Gillis  and  Johnston 
(1983)  for  more  strongly  curved  flows.  An  explanation  for  the 
reversal  in  shear  stress  given  by  Honami  (1980)  was  presented 
in  Gillis  and  Johnston  (1983).  Within  the  curve,  the  dominant 
production  terms  in  the  budget  equation  for  shear  stress  are 
given  by 

For  a  flat  wall,  the  eecond  term  on  the  right  hand  side  is  sera. 
As  the  flow  enters  the  curve,  however,  the  second  term  appears. 
Since  u*5  is  usually  greater  than  o'*,  the  production  rate  is 
decreased  by  this  term,  and  even  goes  negative,  as  shown  by 
calculations  performed  for  the  high  TI  case.  It  is  interesting  to 
note  that  a  large  decrease  in  the  production  level  in  the  study 
by  Gillis  and  Johnston  (1983)  occurred  due  to  the  relatively 
small  radius  of  curvature  ( S/R  =  0.1),  while  a  large  decrease 
in  the  production  level  in  the  present  study  is  obtained  by 
increasing  the  turbulence  intensity  while  keeping  the  radius  of 
curvature  large  (S/R  ■  0.03). 

Profiles  of  PP  are  given  on  Fig.  14.  Values  of  PF  higher 
than  the  corresponding  low  TI  values  are  seen  at  station  1, 
while  an  asymptotic  profile  agreeing  with  the  asymptotic  pro* 
file  for  the  low  TI  case  is  achieved  in  the  curve.  This  dramat¬ 
ically  illustrates  the  dominance  of  curvature  over  turbulence 
intensity.  The  profiles  are  seen  to  agree  well  with  the  wall  heat 
flux  measurements.  Profiles  of  PF  normalised  on  the  extrap¬ 
olated,  effective  value  computed  from  the  wall  heat  flux  are 
shown  on  Fig.  18.  The  values  are  expected  to  approach  unity 
as  the  wall  is  approached. 

The  effect  of  turbulence  intensity  on  t'  and  PF  is  similar 
to  its  effect  on  uV  and  PF.  The  flat-plate  values  are  greater 
than  the  corresponding  low-turbulence  values,  but  this  increase 
is  eliminated  by  the  end  of  the  curve  where  the  valuee  agree 
with  the  asymptotic  profile  of  the  low  TI  case. 

Turbulent  Prandtl  number  profiles  are  shown  on  Fig.  16. 
It  is  interesting  that,  with  high  TI,  values  of  Pr»  become  neg¬ 
ative  in  the  wake  region  of  the  boundary  layer  due  to  the  sign 
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reversal  in  shear  stress  with  no  equivalent  reversal  in  the  sign 
of  the  turbulent  heat  flux,  as  discussed  earlier.  This  strong 
variation  in  Pr,  indicates  that  Reynolds  analogy  is  not  valid. 
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Wetkoda  preaantly  oil  at  for  predicting 
lmlnor-to-ttrbwlant  traneltlon .  too  mtkod  aaploye 
aaplriaally  derived  niiing-langtk-fcypotkeet*  (N-M) 
no  dale  far  tranaltlan  etart.  length  and  potn.  another 
aaploye  aultl -agnation  low ■  Reynol da-nuabar  turbulence 
alaaura  aaeala  morning  that  the  aodale  accurately 
aapture  the  propage  tier-  and  aapllfloatlon  of  asternal 
dleturbaniaa  with  tka  turbulanae  production  torn. 

Thla  aoaona  mtkod  km  ehown  aoae  eucoaoe  raaontly 
ukon  dint urban  aaa  ara  large;  but  nor*  eavalopoant  la 


— dad  before  It  oan  be  put  Into  general  us*. 
Multi-equation  aodels  require  tho  solution  of  ssvsrsl 
diffusion  equations,  e.g.  aonentua,  enthalpy, 
turbulenoe  klnstlo  energy  snd  turbulenoe  dissipation. 
This  son  bo  expansive  snd  ouabersoae.  Prosontly  there 
is  no  cssursnee  that  transition  will  bo  prodletod  noro 
accurately  with  aulti-equatlon  nodols  than  with 
slnplor  nodols.  M.H  nods ling  of  transition,  on  tho 
othor  hand,  is  Inexpensive  and  son  von 1 ant,  but  highly 
onplrlosl.  Bolng  onplrlool,  it  Is  trustworthy  within 
ths  dens  In  of  oxporlnantsl  support,  but  not  outs Ids. 
■ooouns  N.M  nodsling  is  Inexpensive  and  has  oast 
osporinantal  support,  It  has  Its  plaoo.  Ohs  oon elusion 
of  ths  AfOSR/Stanford  Conforonoo  on  Coaplex  Turbulent 
Plows  (Kilns  ot  al. ,  1982)  was  that  no  turbulanoo 
eloouro  nodal  (including  higher-order  nodols)  Is 
oufflolontly  gonoral  and  aoeurato  to  narlt  don In an oo 
over  tho  others. 

Tho  present  paper  Investigates  the  use  of 
ousting  N.R  nodols  to  prod  lot  transition  on 
convex -our ved  boundary  layer  flows  undergoing 
lanlnar-to-turbulon t  transition.  Plrst,  present 
flat-wall  transition  nodols  are  tested  against 
raoantly-takan  data  to  dsteralne  what  oonb Inst  Ion  of 
transition  start,  length  and  path  nodols  appears  to  be 
the  boat  for  this  data  set.  Next,  existing  curvature 
corrections  to  the  N.H  turbulanoo  closure  nodal  are 
tested  against  recent  oon  vox-curved,  fully-turbulent 
boistdary  layer  data  to  deternlne  the  preferred 
correction  for  this  data  set.  Finally,  predictions 
fron  a  oonb 1 nod  tranaltlon-wlth-curvature  nodal  are 
oonpared  to  recent  curved-flow  transition  data.  The 
oonb  1 nod  nodal  has  no  curvature  correction  tern  to  the 
transition  nodal;  therefore,  the  oonparlson  Indicates 
whether  such  a  correction  la  neoessary. 

Objectives  of  the  present  Investigation  are: 

(1)  to  choose  the  best  N.H  flat-wall  transition  nodal, 
based  on  the  data  chosen  for  oonparlson, 

(2)  to  choose  the  best  curvature  oorrectlon  to  the  N.H 
turbulanoo  closure  nodal,  baaed  on  the  data  ohoaen  for 
oonparlson,  and 

(3)  to  test  the  oonblned  nodal  against  the  curved-wall 
transition  data. 

The  pr leery  notlvatlon  of  this  work  la  to  test  present 
onplrlosl  correlations  which  would  be  used  as 
subnodels  for  design  nethods  for  predicting  external 
heat  load  on  the  suction  surface  of  a  highly -oon  toured 
gas  turbine  airfoil  In  the  vicinity  of  boundary  layer 
transition. 

Mixing  Length  Hypothesis  Transition  Models 

HH  nodellng  of  transition  oonslsts  of 
onplrlcslly-bssed  subnodels  for  specifying  transition 
•tart,  transition  region  length  and  the  transition 
path  whloh  Is  usually  expressed  as  an  lnternlttency 
factor,  existing  expressions  for  these  subnodels  are 
given  below. 

Transition  Start  Models 

van  Driest  and  Bluner  (1963)  sssuned  that 
breakdown  of  lanlnsr  flow  occurs  when  the  local 
vortlelty  within  the  bowdary  layer  Is  sufficiently 
Strang  sufficiently  far  fron  the  wall.  The  local 
vortlelty  Reynolds  nunber  la  given  by  Tr  -  (y2/v) 
(dU/dy).  If  the  streanwlse  variation  or  Tr  vs.  y 
profiles  were  non  It  or  ed  in  a  lanlnsr  boisidary  layer, 
one  would  find  that  the  Magnitude  of  the  peak  of  the 
profile  and  the  y-looatlon  of  the  peak  value  increase 
with  the  SQuere  root  of  the  streanwlse  dlstanoe.  When 
transition  begins,  these  two  quantities  begin  to 
la  erases  nore  rapidly,  with  streanwlse  dlstanoe,  than 
In  ths  lanlnsr  boundary  lsyer.  This  observation 
supports  the  van  Driest  and  Rluner  hypothesis  that 


transition  occurs  when  the  local  vortlelty  Reynolds 
nunber  beconea  sufficiently  large  sufficiently  far 
fron  the  wall,  van  Driest  and  Bluner (1963)  next 
applied  the  Pohlhausen  theory  to  find  the  following 
relationship  between  Tr  (at  the  profile  naxlnun), 
Reynolds  nunber  based  upon  boistdary  layer  thickness 
and  the  Instantaneous  Pohlhausen  pressure  gradient 
parsneter,  Tr/Rej  ■1*11.  This  equation  was  then 
used  as  ths  basis  for  an  enplrlcal  fit  of  sxperlnental 
data  with  the  Inherent  aasunptlon  that  transition 
ooeurs  at  a  sufficiently  large  value  of  vortlelty 
Reynolds  nunber,  dependent  upon  pressure  gradient  snd 
free-strean  turbulanoo  intensity.  They  proposed  the 
relationships: 


a.  ffyi  -1  *1.96  ReXj#°-5  T,2  A-0.. 

*x,s 

-  1  -  0.0*85  A  ♦  3.36  Re  T,2  A-  0. 

•«,s 

Brown  and  Nsrtln  (1979)  noted  that  the  van  Driest  and 
Bluner  nodel  predictions  oonpare  well  with  data  when 
the  pressure  gradient  Is  zero  or  adverse  but  results 
were  not  found  to  agree  well  when  there  are  favorable 
pressure  gradients.  This  Is  In  oontrast  to  the 
findings  of  Blair  (19B2)  who  stated  that  the  agreenent 
with  his  zero  and  positive  pressure  gradient  data  was 
excellent. 

Seyb  (1967)  anplrlcally  found  the  following 
expression  for  the  start  of  transition: 


Re 


1000 


e,B  1.2  ♦  70  Ti 


‘“'■o.oioi 


TTt^ 


2.62 


It  Is  based  upon  flat-wall  data  with  -0.09  <  1  <  0.10. 
Rote  that  Reg.s  approaches  zero  as  T;  approaches 
Infinity,  whereas  the  Tollnien-Schlichtlng  Halt  of 
stability  based  upon  linear-stability  theory  suggests 
a  lower  llnlt  of  163.  Brown  and  Martin  (1979)  stated 
that,  of  the  van  Driest  and  Bluner  (1963),  Seyb  (1967) 
and  Hall  and  Clbblngs  (1972)  aodels,  only  the  Seyb 
nodel  ccnpared  favorably  to  data  for  adverse,  zero  and 
favorable  pressure  gradients.  Caugler  (1981)  found 
that  results  for  both  the  Seyb  (1967)  and  Dunhan 
(1972)  aodels  agree  with  turbine  blade  data  of  Lander 
(1969).  Hylton  et  al.  (1983)  selected  the  Seyb  nodel 
as  the  best  oholce  because  It  predicted  transition 
only  on  the  suction  surface  of  their  airfoil,  as  was 
true  according  to  the  data. 

Hall  and  Clbblngs  (1972)  later  correlated 
sxperlnental  neasurenents  to  fors  the  relationship: 
R#a,e  ’  190  ♦  exp  (6.88  -  103  T3).  This  nodel  does 
not  Inolude  the  effect  of  streanwlse  acceleration, 
tdtloh  was  added  In  later  aodels  by  Dunhan  (1972)  and 
Abu-Channan  and  Shaw  (1980);  therefore,  the  Hall  and 
Clbblngs  nodel  was  not  incorporated  Into  the  present 
study. 

Dunhan  (1972)  suggested  the  following  enplrlcal 
expression  for  the  start  of  transition: 


Rs#ji  •  [o.27  ♦  0.73  exp  (-80  T3)](550  ♦  680/D), 

where,  D  •  1  ♦  100  f3  -  211 

This  nodel  Is  an  extension  of  the  Hall  and 
Gibb Inga  (1972)  nodel  with  the  effect  of  pressure 
gradient  added.  Caugler  (1981)  found  favorable 
■greeaent  between  the  results  of  the  Dunhan  nodel  and 
turbine  blade  data  of  Lander (1969). 


Abu-Ohannaa  and  Shaw  (i960)  aodlflad  tha  Hall  and 
Gibb Inga  aodel  ao  that,  aa  free-streaa  turbulanea 
Intensity  Inoraaaaa  to  infinity,  tranaltlon  Reynolds 
nuaber  based  upon  aoaentua  thloknaaa  aayaptotas  to  163. 
They  atatad  that  thla  aodlfloatlon  waa  needed  to  sake 
the  aodel  agree  with  the  Tollalen-Sehllchting  Halt  of 
stability.  They  alao  extended  the  Hall  and  Olbblnga 
(1972)  aodel  to  aooount  for  the  effects  of  pressure 
gradient.  Their  available  data  waa  sore  substantial  , 
than  that  of  Dunhas  (1972)  with  free-streaa  turbulenoe 
lntenaltles  over  the  range  0.1  <  T3  <  5f  and  pressure 
gradients  over  the  range  -0.09  <  1  <  0.06.  The 
resulting  fora  of  their  aodel  1st 

"Vo  • 163  *  *tp  lpu>  '  oToHt  f3l 

where  F(l)  -  6.91  ♦  12.751  ♦  63.6912  when  ISO 

and  F(i)  -  6.91  ♦  2.981  -  12.2712  when  1  >  0 

Priddy  and  Baylcy  (1985)  noted  that  predictions  with 
the  Abu-Ghannaa  and  Shaw  aodel  oaapared  well  with 
their  turbine  airfoil  aeaaureaants  on  the  auction 
surfaoe,  but  that  the  aodel  predicted  a  such  larger 
R«.  ,  value  than  Observed  In  the  experlaent  an  the 
pressure  surfaoe.  They  felt  that  the  discrepancy  on 
the  pressure  surfaoe  was  due  to  the  destabilizing 
oencave  curvature  effect.  It  appears  that  the  prealse 
of  this  paper,  that  no  oonvex  curvature  oorreetion  la 
needed  In  the  start  aodel,  la  born  out  by  the  Priddy 
and  Bay ley  data  set. 

Other  transition  start  aodels  Include  those  of 
Michel  (1951).  Granville  (1953),  Salth  and  Gaaberonl 
(1956),  van  Ingen  (1956).  Crabtree  (1958)  and  Jaffa  at 
al.  (1970).  These  aodels  were  not  lnoluded  In  the 
present  study  partially  because  of  their  age  and 
partially  beoause  of  the  reooaaendatlons  of  Hall  and 
Gibb Inga  (1972),  Brown  and  Martin  (1979),  Gaugler 
(1981).  Blair  (1982)  and  Hylton  at  al.  (1983) 
presented  above.  The  McDonald  and  Fish  (1973)  aodel 
was  not  Included  In  the  present  study  because  It 
requires  solving  the  turbulent  kinetic  energy 
aquation;  therefore,  the  McDonald  and  Pish  aodel  was 
oonsldered  to  not  be  of  the  fora  ehoaen  for  the 
present  study. 

Tranaltlon  Length  Models 

Dhawan  and  Naraalaha  (1958)  showed  that  data  for 
transition  length  can  be  correlated  with; 

Reg  -  5.0  Re,,,0-8,  where  Reg  la  a  Reynolds  nuaber 
baaed  upon  transition  length,  t,  and  ReX|B  is  the 
x-Reynolds  nuaber  of  the  start  of  transition  given  by 
an  appropriate  start  aodel.  The  length,  t,  la  the 
atreaawlse  dlstanoe  between  the  positions  where  the 
Interaltteney  la  0.25  and  0.75.  Consistent  with  this 
la  a  aodlflad  relationship  for  0.0  to  0.99 
Interaltteney  presented  by  Dunhaa  (1972); 

Reg  -  16.8  Re,,.0*8. 

Brown  and  Martin  (1979)  showed  that  the  Dhawan  and 
Naraalaha  aodel  predictions  oaapared  well  with  the 
aaaaureaents  of  Brown  and  Burton  (1978)  over  airfoils 
and  was  better  than  the  Debruge  (1970)  aodel.  ttyltcn 
et  al.  (1983)  seleoted  the  Dhawan  and  Naraalaha  length 
aodel  as  the  best,  by  eoaparlaon  agalnat  their  turbine 
blade  data. 

Debruge  (1970)  presented  the  following  eaplrlcal 
relationship  using  the  fora  of  the  Dhawan  and 
Naraalaha  aodel  along  with  the  expanded  data  base  at 
that  tlaei 


Reg  -  0.0168  Re*,,1*28. 

Chen  and  Thyaon  (1971)  expanded  the  above  fora  to 
include  the  loeal  free-e treaa  Mach  nuaber  finding, 
anpirloallyi 

Reg  -  (60  ♦  *.68  M1*82)  Re,,.0*62. 

Brown  and  Martin  (1979)  suggested  the  following 
eaplrloal  relationships: 

Reg  -  0.91  Re,,.  T,  <  0.001 

Reg  -  0.62  Re,,.  T,  >  0.006 

No  other  K.H  transition  length  aodels  are  known 
by  the  authors. 

Transition  Path  Models 

Dhawwi  and  Naraalaha  (1958)  aeasured  the 
Interaltteney,  the  tlae-fractlon  that  the  flow  is 
turbulant-llke,  during  transition  and  found  a 
wlver sal  function  of  Interaltteney  versus 
dlaenslanleas  transition  length: 

7  -  1  -  exp  (-0.912  n2) 

where  n  la  (x  -  x,)/(x,-  x.)  and  x,  and  x.  are  for 
transition  start  and  end  based  upon  0.25  and  0.75 
Interaltteney,  respectively.  Proa  these  results,  one 
can  write  a  relationship  for  the  Interaltteney  range 
0.0  <  7  <  0.99,  following  Dunhaa  (1972),  aa: 

7  •  1  -  exp(-9.65  n?). 

Thla  path  aodel  was  selected  as  the  best  aodel  by 
Hylton  et  al.  (1983)  arter  eoaparlaon  with  their 
airfoil  data. 

Chen  and  Thyson  (1971 )  drew  rroa  the  theory  of 
Eaaons  (1951)  to  develop  their  transition  path  aodel. 
laaons  round  that  randoa  turbulent  spots  grow  while 
keeping  their  shape  as  they  aove  dowistreaa;  finally 
they  aerge  and  fora  turbulent  boundary  layers.  He 
presented  an  Interaltteney  factor  based  upon  this 
reasoning: 

7  -  1  -  exp(-  gdV*) 

where  g  la  the  souree  rate  density  and  V  is  the 
Influence  voluae  at  a  given  point  for  the  resultant 
turbulent  spot.  The  influence  voluae  at  a  given  tlae 
oan  be  calculated  knowing  the  locations  of  the 
boundaries  of  turbulent  spots  for  tlaes  earlier  than 
the  tine  at  which  the  Interaltteney  la  being  evaluated. 
Chen  and  Thyson  extended  this  aquation  by  evaluating  V 
In  tarns  of  transition  lwigth.  A  fora  of  their 
result,  which  acoounts  for  variations  In  free-streaa 
velocity,  as  given  by  Bradshaw  et  al.  098D  Is: 

7  -  1  -  exp{-G(x  -  x.)  —I 

where  G  la  the  spot  foraatlon  paraaeter.  They  found 
that  the  foraatlon  rate  of  turbulent  spots  depends  not 
only  on  the  transition  Reynolds  nuaber  but  also  on  the 
looal  Mach  nuaber.  With  the  data  available  to  then, 
they  foraulated: 

0-3  u [v^(60  ♦  9.68  ]. 


Mot*  that  the  pressure  gradient  affeets  the  Chen  and 
Thyson  path  aodel  whereaa  the  Abu-Ghannaa  and  Shaw 
(1960)  lndleated  no  pressure  gradient  effect. 

Abu-flhannaa  and  Shaw  (I960)  empirically  found  the 
lateral ttency  relationships 

T  -  1  -  *xp(-5n3). 

their  aeaaureaenta,  supported  by  those  of  Schubauer 
and  Klebanoff  (1959).  showed  that  Interalttanoy  grows 
acre  slowly  for  aasller  distances  beyond  the 
transition  start  and  aore  rapidly  for  largo**  distance*. 
Matching  this  behavior  requires  a  larger  exponent  than 
In  the  Dbawan  and  Baraslaha  aodel. 

■o  other  path  aodels  are  known  to  the  authors. 

Turbulent  flow  Curvature  Modifications 

It  Is  generally  known  that  convex  curvature 
suppresses  turbulence  activity  In  a  turbulent  boundary 
layer  whereas  oanesve  curvature  enhances  activity 
(e.g.  Slaon  and  Moffat  (19»3).  and  lou  at  al.  (1966a 
and  1966b),  So  and  Hellor  (1975)).  The  curvature 
effect  on  turbulent  transport  aust  therefore  be 
Included  into  the  aodellng.  This  effeot  Is  expeoted 
to  be  visible  soon  after  the  beginning  of  transition 
whan  the  Interalttanoy  factor  begins  to  turn  on  the 
turbulent  transport  aodel.  The  curvature  effect  Is 
Introduced  as  s  oor reotlon  to  the  standard  flat-wall 
alxlng  length  turbulence  closure  aodel.  The 
corrections  discussed  herein  take  the  fora  of 
aultlpllera  an  the  alxlng  length  distribution  sueh  as 
the  distribution  given  by  Kays  and  Crawford  (1980). 

Bradshaw  (1973)  aodlfled  the  alxlng  length  for 
the  effeot  of  streaallne  curvature  by  using  an  analogy 
with  ataoapherlo  flows.  He  suggested  that,  for  aaall 
curvature,  the  codification  of  the  alxlng  length  la 
linearly  misted  to  a  particular  fore  of  the 
Richardson  n unbar i 

»o  -  *oO  '  •»!>. 

Where  Rl,  the  Richardson  nuMber,  Is  2U/(R  *U/Jy),  R  la 
the  radius  of  wall  curvature— positive  for  oonvex 
ourvature,  •  Is  an  aaplrloal  oonatant,  te  is  the 
curved-flow  alxlng  length  and  Iq  la  the  flat-wall 
alxlng  length.  Bradshaw  maoaasnded  values  for  •  of 
*.5  and  7.0  for  oonoave  and  oonvex  ourvature, 
respectively.  His  reooaaendatlon  la  based  on  data 
with  alld  streaawlse  ourvature,  6/R  -  0.01.  Later, 
Johnston  and  Elde  (1976)  reoona ended  a  I  value  of  6.0 
for  both  oonoave  and  convex  ourvature  based  on 
ooapressor  blade  conditions.  In  the  present  atudy,  a 
value  of  5.0  was  found  to  be  opt 1 nun,  as  will  be 
dlsouased. 

Adaas  and  Johnston  (1983)  suggested  a 
Modification  to  the  nixing  length  for  oonvex  curvature 
In  both  the  log-linear  region  and  the  outer  region. 
They  found  that  In  the  outer  layer  the  alxlng  length 
noma  Used  on  wall  radius  beoeoes  Invariant  with 
streaawlse  dlstanoc  when  the  ourvature  la  sufficiently 
Strang.  They  east  this  sta tenant  into  the  form 

V*  -  0.0025  tanh  (3*.5  «/*). 

Their  proposed  fore  of  the  near-wall  Modification  for 
strong  ourvature  was  given  as< 

to  •  to/(1  ♦  661) 

where  Rl,  the  Rlohardson  nunber,  1*  given  as  6(1  ♦  8). 
The  stability  parses ter,  6.  la  ooaputed  as  2U/(R  »U/»y 
-  0).  Adana  and  Johnston  rs Banner  dad  a  •  value  of  5.0 
for  the  inner  region  aodel  and  aaployed  a  lag  on  the 


outer  region  aodel  for  situations  where  the  strength 
of  curvature,  «/R,  decreases  abruptly  with  streaawlse 
dlstanoe.  Thla  aodel  laposes  a  first-order  lag  on  the 
outer-layer  modifier! 

where  (te'Oeo  *•  the  local  equlllbrlua  outer-layer 
alxlng  lengths*  given  by  the  above  le/R  equation.  A 
value  of  0.06  was  used  for  the  lag  constant,  C, 
following  their  reoooMendatlon. 

Curvature  Modifications  to  Transition  Models 

There  Is  no  appropriate  ourvature  aodlf lost Ion 
for  M.H  transition  aodellng.  The  three  experlaental 
prograas  dlsouased  below  lend  lnforaatlon  about  this 
rathsr  ooaplex  situation,  however. 

Clauser  aid  Clauaer  (1937)  aeasumd  transition 
looatlona  on  oonoave  and  oonvex  walls.  They  concluded 

that  Uta  transition  start  x -Reynolds  nunber  ohanges 
linearly  with  xs/R.  Proa  their  data  a  mlatlonshlp  can 
; be  foraulated  ass 

■*x,s  “  (  *.50  ♦  0.25  x,/R  )  *105. 

The  relationship  gives  a  delay  of  transition  on  the 
oonvex  wall  and  an  early  transition  on  the  oonoave 
wall.  The  free-streaa  turbulence  level,  now  known  to 
be  an  laportant  par ana  ter,  was  not  given;  but,  fron 
their  data  and  riat-wall  lnforaatlon  now  available, 
one  would  deduoe  that  their  turbulanoe  Intensity  was 
around  on*  peroent.  On*  orltlclaa  of  their 
oonolualon  Is  the  near  linearity  which  they 
olaln— one  would  expect  the  transition  ae chan lsa  on 
the  oonoave  wall  to  be  fun daaen tally  different  than 
that  on  the  oonvex  wall  and,  therefore,  the  effect  on 
transition  start  would  be  different.  Also,  their 
experiaent  does  not  lnoorporat*  the  cross-effects  of 
ourvature,  free-streaa  turbulanoe  Intensity  and 
streaawlse  acceleration.  The  cross-effect  of 
ourvature  and  free-streaa  turbulanoe  intensity  was 
Shown  by  Hang  and  81aon  (1965)  to  be  significant. 

Llepaann  (19*3)  Measured  transition  start 
locations  on  oonoave  and  oonvex  walls.  The  documented 
free-streaa  turbulanoe  level  was  0.1  percent.  He 
found  a  strong  effect  of  oonoave  curvature  and  no 
effect  of  oonvex  ourvature  on  transition.  The  finding 
that  oonvex  ourvature  has  no  influence  on  the 
transition  start  location  for  lower  turbulence 
Intensity  runs  is  oontrary  to  the  findings  of  Clauser 
■id  Clauser  (1937)  and  Hang  and  Slaon  (1965). 

Llepaann  notes  that  in  his  curved-flow  runs, 
turbulanoe  froa  the  opposite  (oonoave)  wall  la  spread 
to  the  oonvex  wall  raising  the  "free-streaa" 
turbulanoe  intensity  value  froa  0.1  peroent  at  the 
beginning  of  the  test  plat*  to  0.75  peroent  at  the  end. 
This  Increased  turbulanoe  Intensity  over  that  of  the 
Llepaann  straight-wall  measurements  would  result  In  an 
earlier  transition  on  both  oonvex  and  oonoave  walls. 
This  earllar  transition,  superposed  with  a  retarding 
effeot  of  oonvex  curvature  and  an  enhancing  effect  of 
oonoave  ourvature  nay  have  led  to  an  inoorrect 
oonolualon  about  the  ourvature  effeot— the  ooncluslon 
that  oonvex  curvature  has  no  effect  on  the  transition 
start  location  and  oonoave  curvature  has  a  large 
effect. 

Hang  and  Slaon  (1965)  Bad*  transition  boundary 
layer  host  transfer  and  fluid  Mechanics  Measurements 
an  a  wall  that  was  straight,  then  bent  to  two 
different  radii  of  convex  curvature.  Tests  were  run 
at  0.66  poroant  and  2  peroent  free-streaa  turbulanoe 
Intensity.  They  foisid  a  doubling  of  the  transition 


■tart  location  fro*  tha  fiat  wall  eaaa  to  a  eiidly 
curved  aaaa  whan  tha  turbuleno*  Intensity  waa  0.68 
paroant ,  but  no  furthar  lnoraaaa  in  tha  affaet  of 
eurvatur*  aa  tha  curvature  strength  waa  Increased 
furthar.  Thay  alao  found  only  a  a  In  or  af fact  of 
aurvatura  on  tha  tranaltlon  atart  position  at  tha 
Iklghar  fraa-atraaa  turbulanoa  level. 

laoauaa  af  tha  atreng  aroaa ■ aerralatlen  of 
fraa-atraaa  turbulanoa  affaota  and  aurvatura  af faota , 
aa  doeu— ntad  by  Nang  and  tlaon  (IMS),  and  baoauaa  of 
tha  aonfllating  aonolualona  froa  tha  above  thraa 
expert— ta.  It  la  dlffiault  to  aenaolldoto  tha  abora  ( 
rooulta  Into  a  curvature  aorroatlon  far  tanaltlan 
nodal  a.  Tha  hlghar  turbulanoa  aaaa  of  Nang  and  (loon 
auggaata  that  for  hlghar  turbulanoa  larala,  a 
aurratura  oorraotlon  to  tha  tranaltlon  nodal  nay  not 
ba  naoaaaary,  howarar.  Thla  la  a  hypothaala  whlah  tha 
proa ant  papar  taata. 


SCLCTIOB  MOOBUB 

Tha  tranaltlon  atart,  langth  and  goth  nodal*  and 
tha  turbulent  aurratura  nodlfloatlana  dlaauaaad  abora 
wora  lnoorporatad  Into  tha  boundary  layer  coda  STANS 
(Crawford  and  Kays,  1976).  Laminar  raloolty  and  total 
enthalpy  profile*  war*  given  at  the  atart in g  looatlon 
1  an  downatraan  of  tha  loading  adg*  of  tha  taat  plata. 
for  oaaaa  which  war*  Initiated  aa  turbulent  boundary 
layers,  naaaured  velocity  and  total  enthalpy  profile* 
war*  auppllod  for  the  beginning  of  oonputatlon.  Tha 
naaaured  wall  aurfaoa  boat  flux  variation  with 
atroanwla*  dlatanoe  was  used  as  the  thermal  boundary 
condition.  Moat  cases  ware  run  aaaunlng  uniform 
properties  across  the  boundary  layer  after  test  cases 
war*  run  to  verify  that  variable  property  effeota  were 
Insignificant.  The  turbulanoa  closure  aodel  used  waa 
the  nixing  langth  nodel  onploylng  a  standard  nixing 
length  distribution  as  dlacuaaed  In  Kays  and  Crawford 
(1980).  Turbulence  oloaure  In  the  energy  aquation  was 
by  use  of  specified  turbulent  Prandtl  n unbar a.  Those 
turbulent  Prandtl  nun  bars  war*  the  vaaaured  values, 
except  In  the  eases  of  Blair  and  Merle  0  981,  1982) 
where  a  oonatant  turbulent  Prandtl  nuaber  of  0.90  was 


The  transition  start  nodals  used  different 
fraa-atraaa  turbulanoa  intensity  definitions,  Tj  and 
T3.  When  available,  T*  was  used:  when  necessary,  T, 
waa  substituted.  The  local  free-streaa  turbulence 
lntanaity  wa*  linearly  interpolated  froa  the  data 
points. 


BKSMLTS  AMD  DISCUSSION 

Bel act Ion  of  Test  Casas 

Thro*  sxparlnantal  studies  (Vang  at  al.,  1965, 

You  at  al. ,1986a  and  1986b  and  Vang  and  Slacn,  1985), 
all  don*  at  tha  Neat  Transfer  Laboratory  of  the 
University  of  Minnesota,  war*  chosen  aa  the  primary 
ooaparlaon  oases.  These  three  war*  den*  on  the  saa* 
facility  using  the  aaaa  experimental  procedures. 

Thay,  therefor*,  represent  a  coherent  data  aet  where 
aurvatura  and  transition  af foots  are  first  Isolated, 
than  combined.  Those  tests  are  familiar  to  the 
authors  and  data llad  data  froa  tha  teats  are  available 
to  tha  authopi. 

•lx  eases  of  81alr  and  Varle  (1981.1982)  were 
alao  ahoaan  to  lnveatlgate  the  affects  of  fraa-stroam 
turbulanoa  intensity  and  atroaawisa  pressure  gradient 
on  tha  transition  atart. 


The  Tost  Cases 

An  experimental  study  of  transition  on  a  flat 
wall  by  Vang  at  al.  (1985)  was  ehosen  a*  a  test  oas* 
for  dateralnlng  the  boat  transition  start,  langth  and 
path  submodels.  In  their  study,  data  for  0.66  paroant 
and  2  paroant  turbulanoa  intensities  with  xero 
atraamwlsa  proaaur*  radiant  war*  praam* ted.  Tha  wall 
beat  flux  waa  nearly  tnifera.  The  2  paroant 
turbulanoa  Intensity  ease  was  aaleetad  aa  tha  primary 
aaaa  for  comparison.  It  waa  shown  to  b* 
two-dimensional  throughout  tha  taat  langth  and  good 
momentum  and  anargy  balances  on  tha  taat  wall  war* 
achieved.  Tha  measured  turbulanoa  intensity  varied 
from  2.3  paroant  at  the  loading  adg*  to  1.9  paroant  at 
the  downstream  and.  Though  the  akin  friction  was 
measured  at  only  nine  atraamwlsa  locations,  local  host 
transfer  coefficient  data  was  taken  at  more  than  50 
atraamwlsa  positions  making  tha  transition  path 
clearly  dlsoemabl*.  Comparison  runs  wars  alao  mad* 
ualng  tha  0.68  paroant  turbulanoa  Intensity  oasa  of 
Hang  at  al.  (1965)  and  with  tha  data  of  Blair  and 
Marla  (iMl,  1982)  with  and  without  favorable  pressure 
grad 1 ants  and  for  several  values  of  fraa-atraaa 
turbulanoa  Intensity  using  tha  submodels  aalacted. 

The  aurvatura  oorraotlon  to  the  turbulent  mixing 
langth  was  chosen  by  oamparlaan  with  the  experimental 
data  of  You  at  al.  (1986a,  1966b)  measured  in  a  fully 
turbulent  flow  of  moderate  oonvex  curvature  strength 
(0.013  <  t/K  <  0.03).  In  this  experiment  a  turbulent 
bowdary  layer  was  grown  on  a  flat  wall  followed  by  a 
constant-radius  curved  section  and  a  downstream  flat 
wall.  The  data,  therefore,  tests  tha  model's  response 
to  both  the  Introduction  of  and  recovery  from  oonvex 
eurvatur*.  Two  values  Tor  the  strength  or  curvature 
war*  achieved  by  banding  the  test  wall  to  two 
different  radii.  Comparisons  were  mad*  with  eases  of 
0.65  peroent  and  1.85  paroant  turbulence  intensity. 

The  combined  model,  with  the  ehosen  transition 
start,  langth  and  path  submodels  aa  well  as  the 
eurvatur*  modification  to  the  turbulanoa  nodel,  was 
tasted  against  the  transitional  flow  data  of  Vang  and 
Simon  (1965)  taksn  on  a  oonvex  wall.  Comparisons  were 
mad*  for  0.68  peroent  and  2  paroant  turbulence 
Intensities,  each  at  two  radii  of  streanwis* 
curvature. 

Tranaltlon  on  flat  Vails 

A  comparison  of  predicted  stremswls*  variations 
of  Stanton  number  on  a  flat  wall  using  the  transition 
atart  models  discussed  earlier  la  shown  on  figure  1 (m) . 
Bach  run  was  made  with  the  same  transition  length  and 
path  models.  Shown  also  is  the  data  from  the 
comparison  experiment  discussed  above.  The 
Abu-Ghannam  and  Shaw  (i960)  atart  model  accurately 
predicted  tha  looatlon  of  tha  onset  of  transition 
whereas  the  other  three  models  gsve  s  later  transition 
than  the  data  lndleata.  A  oamparlaan  of  further  test 
results  of  the  start  models  against  various  degrees  of 
frea-stream  turbulence  lntanaity  and  streamwls* 
pressure  gradient  Is  shown  In  Table  1 .  The  Abu-Ghannam 
and  Shaw  (i960)  start  model  appears  to  be  better  than 
tha  other  atart  models.  Based  upon  those  comparisons, 
and  favorable  results  of  subsequent  comparison*  using 
the  Abu-Ghannam  and  Shaw  atart  model  and  the  ehosen 
length  and  path  models,  the  Abu-Ghannam  and  Shaw  start 
aodel  was  ohosen. 
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Flguro  1(b)  shows  tha  aaoa  data  and  pradlotlons 
using  ths  AbwOhannan  aid  Shaw  start  nodal  and  tha 
various  langth  nodala  dlaouaaad  aarllar.  All  thraa 
pradlotlons  onploy  tha  Sana  transition  path  nodal. 
All  thraa  nodala  pradlet  s  ahortar  transition  langth 
than  tha  data  lndloat*.  Tha  Brown  and  Martin  (1979) 
nodal  gava  alaost  tha  sans  rasults  as  tha  Dabrugo 
0970)  nodal.  Tha  Dhawan  and  Maraalaha  0958)  nodal 
glass  tha  bast  raaulta ,  however.  Baaad  an  this 
aangorlsan  and  favorable  raaulta  of  subsequent 
aanparlsons  using  tha  AbvOhaman  and  Bhaw  start 
nodal,  tha  Dhawan  and  Baraaltfta  langth  nodal  and  tha 
ahaoan  path  nodal,  tha  Dhawan  and  Raraslaha  langth 
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Raaulta  of  pradlotlons  of  tha  aant  data  sat  using 
tha  thraa  transition  path  nodala  ars  shown  on  fig. 
He).  Tha  oonparlson  ossa*  wars  all  oonputad  with  tha 
AbwQhannan  and  Shaw  start  nodal  and  tha  Dhawan  and 
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lirulihi  length  aodel.  The  Chan  and  Thyaon  (1971) 
aodel  la  elaarly  passing  to  turbulent-llke  flow  too 
qulokly.  Tha  Dhawan  and  II araslaha  (1958)  and 
Abu-Ghannaa  and  Shaw  (i960)  aodala  ahow  soaewhat  tha 
aaaa  path  but  pradlotad  Stanton  nuabers  froa  tha 
Abu-Ghannaa  and  Shaw  aodal  ara  alowar  to  ohanga  with 
downatraaa  dlatanoa  than  thoaa  of  tha  Dhawan  and 
Haraslaha  aodal.  Thla  alowar  aarly  rasponsa  of  tha 
Abu-Ghannaa  and  Shaw  aodal  la  saan  in  tha 
lntaralttanoy  Tarsus  dlstanoa  funotlons  shown  on  Fig. 
2.  Baeausa  of  thla  aarllar  rlaa  In  lntaralttanoy  and 
baeausa  of  favorable  subsequent  ooaparlaona  of  tha 
raaulta  of  tha  ooablnad  aodal  to  experimental  data, 
tha  Dhawan  and  H araslaha  (1956)  aodal  was  ehoaan. 

Tha  transition  start,  langth  and  path  aodala  hava 
now  bsan  aalaotad.  Coaputad  flat-wall,  transition  akin 
frlotlon  and  Stanton  nuabar  raaulta  ara  ooaparad  to 
tha  data  of  tha  ooaparlaon  tast  on  Pig.  3.  In  Pig.  A, 
raaulta  froa  tha  ehoaan  ooablnatlon  of  transition 
aubaodals  ara  ooaparad  to  1  percent,  2  percent ,  and  A 
percent  turbulence  intensity  data,  respectively ,  froa 
Blair  and  Karla  (1962).  Pig.  5  shows  predictions,  by 
tha  sase  aodal,  of  the  Blair  and  Kerle  (1961)  data  for 
ooablnatlona  of  free-strean  turbulence  Intensity  and 
streaawlse  favorable  pressure  gradient.  A  favorable 
pressure  gradient,  dp/dx  <  0,  stabilises  the  flow 
delaying  transition!  an  adverse  pressure  gradient  has 
the  opposite  effect.  Of  all  the  botatdary  layer  fluid, 
tha  low  aoaantua  near-well  flow  tends  to  be  aost 
affected  by  pressure  gradient.  Thla  la  also  the 
region  of  aaxlaua  vortlelty  Reynolds  nuabar  and  the 
region  aoat  likely  to  aapllfy  aaall  disturbances 
eventually  leading  to  transition.  A  ooaparlson  of  the 
data  of  Pig.  5  shows  the  stabilising  effect  of  a 
favorabla  pressure  gradient.  Sensitivity  to  pressure 
gradient  would  be  greatly  reduced  if  the  near-wall 
flow  possessed  aore  streaawlse  aoaentua.  This  oan  be 
affected  by  Increasing  the  free-streaa  turbulence 
Intensity.  A  ooaparlson  of  Pigs.  A(e)  and  5(o)  ahows 
a  aueh  reduced  effect  of  aooeleratlcn  froa  that  of  the 
previous  oeaparlson  (Pigs.  5(b)  and  A(b)).  In  the 
experlaant  of  Abu-Ghannaa  and  Shaw  (i960),  the  effect 
of  pressure  gradient  nearly  disappeared  for  turbulanoe 
Intensities  larger  than  5  percent.  The  pressure 
gradient  and  turbulanoe  Intensity  effaots  of  the  Blair 
>nd  Karla  (1961 )  data  are  captured  reasonably  well  by 
the  chosen  aodal. 


Fig  4  Transition  with  zero  pressure  gradient 
(•at  piste) 
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Figure  6  shows  a  oonparlson  of  predictions  and 
the  Stanton  nunber  data  of  the  flat-wall,  lower 
turbulence  Intensity  (0.66f)  oase  of  Wang  et  al. 
(1985).  The  predloted  onset  o?  transition  agrees  well 
with  the  neasurenents  but  a  large  discrepancy  la  seen 
between  the  transition  path  data  and  predictions, 
tang,  et  al.  noted  that  the  naaantus  and  energy 
balanoee  used  to  verify  the  test  began  to  degrade  at 
the  start  of  transition  due  to  the  development  of 
three-dlnanslenallty  caused  by  earlier  transition  on 
the  aide-walls.  This  is  thought  to  not  account  for 
the  full  discrepancy  between  the  data  and  the 
predictions,  however.  To  the  authors'  knowledge,  no 
detailed  transition  path  data  In  addition  to  the  Vang 
et  al.  (1985)  oase  are  available  In  this  turbulence 
Intensity  range. 


Fig.  7  Comparison  oi  turbulent  curvature  modificationt 

(a)  T,  •  0  65  V  R  •  SO  cm  (oortvei) 

(b)  T,  -  0.65  %,  R  •  210  cm  (convex) 

(c)  R.SOem(oonvei) 

Turbulent  Flow  on  Convex-Curved  Walls 

The  two  eurvature  nodlflcation  nodels  previously 
discussed  were  tested  against  the  experlnental  data  of 
You  et  al.  (1986a,  1986b).  The  results  are  shown  In 
Fig.  7.  Both  nodels  predict  the  three  oases 
reasonably  well  but  the  Adana  and  Johnston  (1983) 
nodal  is  sonewhat  better,  particularly  In  its 
perfornanoe  against  the  rsoovery  data  on  the  flat  wall 
downstrean  of  the  curve.  The  Adana  and  Johnston 
eurvature  nodlfler  was  therefore  chosen. 

Transition  on  Convex-Curved  Valla 

The  eonblned  nodel,  consisting  of  the  Abu-Channan 
and  Shaw  (1961)  start  nodel,  the  Dhawan  and  Naraslnha 
(1956)  langth  and  path  nodels  and  the  Adans  and 
Johnston  (1963)  eurvature  nodlfler  was  next  tested 
against  curved-wall  transition  data.  The  oonparlson 
oases  were  the  experlnants  of  Vang  and  Slnon  (1965) 
taken  with  two  radii  of  eurvature,  each  with  two 
levels  of  turbulence  intensity  (0.66  and  Zl).  The 
oon par Isons  with  the  two  percent  turbulenoe  eases. 
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shown  on  rig.  8,  are  quit*  favorable.  Under  these 
eondltlona,  It  appear*  that  there  need  not  be  a 
curvature  oorreotlon  to  the  transition  start 
nodel— the  curvature  Modifier  to  the  turbulence  aodel, 
uhlah  would  alter  the  transitional  and  turbulent  flow 
data  but  not  the  transition  start  location,  aeons  to 
aooount  for  the  entire  ourvatur*  of foot.  The 
oaaparlscn  against  the  low  turbulsno*  Intensity  data 
(Fig.  9)  is  not  so  favorable,  however.  Clearly,  for 
low  turbulenoa  intensity,  a  different  transition 
nodel,  which  Includes  a  ourvatur*  oorreotlon,  la 
needed.  Though  the  turbulenoa  Intensity  value  at 
whloh  suoh  oorreotlon  to  the  aodel  Is  needed  Is  not 
preolsely  known.  It  would  appear  to  be  sonewher* 
around  on*  to  two  peroent.  The  ooncluslon,  then,  la 
that  the  nodel  chosen  in  the  present  study  can  be  used 
to  predict  transition  an  ourved  walls,  provided  that 
the  free-strean  turbuleno*  intensity  value  Is  about 
two  peroent  or  greater.  This  Is  useful  because  gas 
turbine  suet Ion  surface  design  eondltlona  would 
•morally  not  violate  this  requlranant.  I 


ccgausigs  and  gcowgMMMOgs 

The  following  conclusions  ean  be  drawn  fran  the 
present  study: 

1.  A  transition  nodel  based  on  the  Abu-Ghannaa 
and  Shaw  transition  start  subnodel,  and  the 
Dhawan  and  Haras Inha  length  and  path 
subnodela  la  the  best  transition  aodel  of  thr 
nod* Is  tested  above  for  oonputing  the 
behavior  of  the  chosen  flat-wall,  transition 
oonparlson  case  data.  Predictions  with  this 
nodel  were  shown  to  also  be  favorable  when 
oonpared  to  transition  data  that  Included 
various  levels  of  free-strean  turbuleno*  and 
favorable  streanwls*  pressure  gradient 
especially  when  the  free-strean  turbuleno* 
intensity  value  renained  higher  than  about 
an*  peroent. 

2.  The  ourvatur*  Modification  to  the  nixing 
length  nodel  proposed  by  A  dans  and  Johnaton 
was  the  preferred  aodel  for  predicting  the 
effeot  of  oonvex  ourvatur*  as  displayed  In 
the  ohoaen  turbulent  ourved  flow  oase. 

3.  The  oonblned  (transition  and  ourvatur*)  aodel 
predicts  the  eonvexly  ourved  wall  transition 
data  well  for  the  high  turbuleno*  Intensity 
oase  (21) ,  but  a  ourvatur*  oorreotlon  to  the 
transition  nodel  la  needed  for  the  low 
turbuleno*  oase  (0.681). 

Jn  general,  nor*  oonparlson  data  fron  ourved 
transitional  flows  would  be  useful  for  oohtlnued  aodel 
developnent.  Nor*  renalna  to  be  don*  on  predicting 
oonvexly  ourved  transitional  flows  with  various 
degrees  of  turbuleno*  Intensity.  Data  la  needed  to 
support  the  present  paper's  prellalnary  finding  that 
no  ourvatur*  oorreotlon  to  the  transition  aodel  is 
needed  when  the  free-streaa  turbuleno*  Intensity  Is 
high  (22S),  and  for  Modeling  transitional  flows  with 
low  turbuleno*  intensity.  The  effeot  of  oonoav* 
curvature  on  transition  la  presently  not  well 
understood— sore,  detailed  transition  path  data  In 
this  very  ooaplex  flow  is  needed. 
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ABSTRACT 

Measurements  of  film  cooling  on  a  simulated  turbine  blade  are 
conducted  using  a  mass  transfer  technique.  Under  the  influence  of 
the  endwall,  dramatic  changes  of  film  cooling  performance  occur  on 
the  convex  surface  of  the  blade  as  compared  to  the  region  where 
the  flow  Is  two-dimensional.  The  result  is  a  triangular  region, 
where  coolant  is  swept  away  from  the  surface  by  the  passage  vortex 
present  between  adjacent  blades.  In  order  to  predict  the  area  of 
this  unprotected  region,  the  Influence  of  several  parameters 
Including  density  ratio,  blowing  rate  and  number  of  rows  of 
Injection  .ioles  are  studied.  The  presence  of  the  endwall  affects 
the  film  cooling  performance  on  the  concave  surface  only  slightly. 


Introduction 


To  Improve  the  performance  of  a  gas  turbine  engine,  the  temperature  of 
combustion  gas  entering  the  first  stage  turbine  has  Increased  to  as  high  as 
1750  K  which  Is  above  the  allowable  limit  for  present  alloys,  used  as  turbine 
materials.  An  effective  cooling  system  Is  required  to  protect  the  gas  turbine 
blade  from  high -temperature  failure.  Surface  film  cooling  combined  with 
Internal  convection  cooling  and  Impingement  cooling  Is  widely  used  In  the 
aircraft  Industry.  A  better  understanding  of  this  cooling  system  can  Improve 
the  design  of  gas  turbine  engines.  The  present  work  focuses  on  film  cooling 
using  a  plane  simulated  turbine  blade  cascade.  Particular  attention  Is 
directed  to  the  Influence  of  the  endwall  on  the  adiabatic  wall  effectiveness. 

Adiabatic  wall  effectiveness  Is  defined  by: 

naw  «  TgLlIr  (1) 

T2  -  'r 

Based  on  the  definition  of  the  adiabatic  wall  temperature,  the  heat  transfer 
Is  defined  by  the  equation: 

qw  «  h  (  Tw  -  Taw)  (2) 

Values  of  the  heat  transfer  coefficient  are  generally  close  to  those  without 
coolant  Injection  If  the  region  near  the  injection  holes  is  excluded. 

The  present  study  applies  the  heat-mass  transfer  analogy  to  Investigate 
the  film  cooling  effectiveness.  Detailed  discussion  about  heat-mass  transfer 
analogy  Is  available  in  ref[l].  A  gas  mixture  containing  a  tracer  Is  injected 
Into  the  mainstream.  The  concentrations  of  the  tracer  are  measured  In  samples 
drawn  through  taps  which  are  located  on  the  surface  of  a  blade.  The  local 
impermeable  wall  effectiveness,  ^Iw,  Is  defined  by: 
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As  the  concentration  of  the  tracer  In  the  Mlnstreaai  Is  zero,  the  Impermeable 
nail  effectiveness  can  be  calculated  fro*  the  equation: 


Early  studies  [2,  3]  of  the  fll*  cooling  on  a  turbine  blade  Mere 
conducted  In  the  central  space  of  the  blade  to  *1n1«1ze  the  Influence  of  the 
enduall.  The  effect  of  b1aden«a11  curvature  to  the  performance  of  the  fll* 
cooling  Mas  reported.  However,  the  span  of  *ost  blades  Is  not  long  coopered 
to  the  chord  length.  The  flow  field  as  Mali  as  the  fll*  cooling  are  strongly 
affected  by  the  presence  of  the  endMall.  A  recent  Mork  [A]  found  that  there 
Is  a  close  relation  betMeen  the  floa  field  In  the  endMall  region  and  the  fll* 
cooling.  Measurenents  of  the  fll*  cooling  In  a  plane  cascade  shorn  that  a 
passage  vortex  travelling  fro*  the  adjacent  blade  smeeps  the  coolant  fro*  the 
convex  surface  and  produces  a  trlagular  region  Mlthout  fll*  cooling 
protection.  An  understanding  of  the  floM  field  In  the  near  enduall  region  Is 
essential  when  conducting  a  study  of  fll*  cooling.  Studies  of  the 
three-dimensional  flow  field  Mlthln  a  turbine  cascade  passage  have  been 
described  by  others  [5,  6,  7,  8].  A  review  of  Investigations  of  the  flow 
field  for  past  decade  Mas  reported  In  ref.  [9]. 

The  present  work  is  an  extension  of  an  earlier  study  done  with  a  single 
row  of  holes,  a  single  blowing  rate  and  a  single  density  ratio.  Here  the 
effect  of  density  ratio,  two  rows  of  Injection  holes  and  blowing  rate  on  film 
cooling  near  an  endwall  are  considered. 


Experimental  System  and  Test  Conditions 


Measurements  arc  conducted  In  an  open  cycle,  low  speed  wind  tunnel  using 
a  plane  cascade  of  six  turbine  blades  [1].  Two  central  blades  with  Injection 
boles  face  each  other.  Four  others  are  solid.  Figure  1(a)  shows  the  location 
of  the  Injection  holes  and  sampling  taps.  Two  staggered  rows  of  Injection 
holes  Inclined  at  an  angle  of  35*  are  used.  The  angles  of  flow  Inlet  and 
outlet  with  respect  to  the  blades  are  shown  In  Fig.  1(b).  To  assure  uniform 
flew  though  each  hole,  the  Injection  hole  half  embedded  in  the  endwall  Is 
plugged.  The  lateral  positions  of  each  of  the  four  rows  of  sampling  taps  are 
2/0  •  0.0,  0.5,  1.0,  and  1.5,  respectively.  Each  row  has  seven  sampling  taps 
distributed  at  various  positions  downstream.  The  central  two  blades  can  be 
moved  In  or  out  of  the  endkall  using  a  positioning  mechanism  underneath  the 
blades.  For  any  single  run,  the  blades  are  positioned  at  a  fixed  distance 
from  the  enfeall,  that  Is,  the  lowest  hole  of  the  second  row  Is  centered  on 
the  en*»al1.  For  each  position,  samples  are  drawn  from  all  four  rows  of 
sampling  taps.  These  represent  the  regions  between  different  paired  holes, 
that  Is,  only  one  row  of  sampling  taps  lie  Inside  the  shaded  area  where  the 
meesurement  region  Is.  To  get  data  for  results  between  a  given  pair  of  holes 
(at  fixed  H/0),  four  different  posltlonlngs  of  the  blades  are  required.  For  a 
given  blowing  rate  and  a  density  ratio,  the  blades  will  be  positioned  Into  a 
dozen  different  locations  to  conduct  a  series  of  measurements  before  the 
two-dimensional  flow  region  Is  reached.  Note  that  the  shaded  area  only  covers 
half  the  space  between  a  given  pair  of  holes  as  the  sampling  taps  are  designed 
for  a  study  In  the  central  span  of  the  blades  where  the  flow  Is  symmetrical 
around  the  Injection  hole. 

The  Injected  flow  Is  produced  by  mixing  refrigerant -12  or  helium  and  air. 
Two  density  ratios,  R  ■  2.0  and  0.96,  are  used.  For  the  measurments  on  a 
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convex  surface,  three  blowing  rates,  M  ■  0.5,  1.0,  and  1.5,  are  used  for  both 
density  ratios.  On  the  concave  surface,  the  blowing  rates  are  1.38  and  1.5 
corresponding  to  the  density  ratio  of  0.96  and  2.0,  respectively. 

The  sampled  gases  are  analyzed  using  a  gas  chromatograph.  The  velocity 
profile  Is  measured  14  cm  upstream  of  the  leading  edge  of  the  blades.  The 
free  stream  velocity  Is  10.2  */s.  The  boundary  layer  Is  turbulent  with  a 

thickness,  0f  m  corrt$p0n<n ng  to  H/g  ■  4.2.  The  turbulent  Intensity 
of  the  free  stream  Is  1.2%.  The  momentum  thickness  Is  1.2  mm  and  the 
displacement  thickness  Is  1.8  mm.  By  calculating  a  virtual  origin  of  the 
turbulent  boundary  layer  and  assuming  a  variation  of  a  turbulent  boundary 
layer  on  a  flat  plate,  the  boundary  layer  thickness  at  the  location  of  the 
stagnation  line  of  the  blade  would  be  about  6%  larger  than  the  value  measured 
14  cm  upstream  of  the  leading  edge  of  the  blade.  The  mainstream  velocity  near 
the  Injection  holes  of  the  first  row  Is  17.8  ”/s  on  the  convex  surface  and  4.2 
m/%  on  the  concave  surface.  The  uncertainty  of  measuring  Impermeable  wall 
effectiveness  Is  6%  as  the  value  of  the  effectiveness  Is  greater  than  10%. 

The  uncertainty  Is  7%  In  calculating  the  blowing  rate  and  4%  In  calculating 
the  density  ratio. 

Flow  Field 

Because  of  a  strong  relation  between  film  cooling  and  the  flow  field,  a 
detailed  description  of  the  flow  field  Is  essential  to  explain  the  results  of 
the  film  cooling  measurements.  A  sketch  of  the  flow  field,  based  on  results 
of  flow  visualization  [3],  Is  shown  In  Fig.  2.  The  main  structure  Is  clear 
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but,  without  more  Investigations,  some  details,  especially  close  to  the  blade, 
are  still  ambiguous.  In  Fig.  2,  the  attachment  (stagnation)  line.  Indicated 
by  ai  -  a2,  extends  from  the  stagnation  point  on  the  blade  Into  the 
mainstream.  The  separation  line.  Indicated  by  si  -  s2,  lies  In  front  of  the 
leading  edge  of  the  blade.  A  saddle  point,  Ssi,  Is  a  singular  separation 
point  and  Is  located  at  the  Intersection  of  the  attachment  line  and  the 
separation  line.  Several  flow  regions  are  bounded  by  these  two  separation 
lines  and  two  attachment  lines.  When  the  mainstream  approaches  the  leading 
edge  of  the  blade,  the  Inlet  boundary  layer  separates  and  a  horseshoe  vortex 
Is  generated  behind  the  separation  line.  After  the  mainstream  travels  Into 
the  blade  passage,  a  pressure  variation  Is  produced,  high  pressure  near  the 
concave  surface  and  low  pressure  near  tne  convex  surface.  Behind  the 
separation  line,  one  leg  of  the  horseshoe  vortex  on  the  concave  side 
travelling  downstream  Is  driven  by  the  pressure  difference  In  the  blade 
passage  and  becomes  part  of  the  passage  vortex.  The  other  leg  of  the 
horseshoe  vortex,  on  the  convex  side.  Is  confined  near  the  blade  wall  and 
rolls  up  when  It  encounters  the  passage  vortex  from  the  the  adjacent  blade. 

On  the  endwall,  a  boundary  layer  Is  driven  by  the  pressure  difference  from  the 
concave  side  to  the  convex  side. 

Essentially,  the  passage  vortex  Is  formed  from  the  concave  side  leg  of 
the  horeshoe  vortex  and  fed  from  the  cross  flow  on  the  endwall  and  the 
mainstream  In  the  blade  passage.  It  grows  bigger  as  It  travels  downstream. 
When  the  passage  vortex  encounters  the  convex  surface  of  the  adjacent  blade, 

It  apparently  lifts  up  the  convex  side  leg  of  the  horseshoe  vortex  from  that 
blade.  Also,  a  counter  vortex  Is  generated  underneath  the  passage  vortex. 
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Near  the  endtall,  the  Uniting  streamlines  near  the  blade  surface  are  not 
two-dimensional.  On  the  concave  surface,  the  Uniting  streamlines  are  skewed 
downward  because  fluid  feeds  Into  the  crossflow  on  the  endwall.  On  the  convex 
surface,  the  limiting  streamlines  are  pushed  up  by  the  passage  vortex. 

Initial  work  on  film  cooling  near  the  endwall  [4]  indicates  that  the  passage 
vortex  sweeps  the  coolant  flow  from  the  convex  surface  and  there  Is  a 
triangular  area  In  which  the  coolant  Is  essential  absent. 

Impermeable  Mall  Effectiveness 

Convex  surface 

For  M  •  1.0  and  R  ■  2.0,  Iso-effectiveness  lines  of  the  convex  surface 
are  shown  In  Fig.  3.  These  lines  are  constructed  In  the  measurement  regions 
which  start  at  x/d  •  1*88  and  end  at  at  N/p  •  70.93.  Note  that  these 
measurement  regions  only  cover  half  area  between  two  adjacent  Injection  holes. 
For  the  region  of  31.5  >  H/p  >  30,  the  coolant  jet  is  symmetric  around  the 
Injection  hole  centerline  as  can  be  seen  from  the  periodicity  of  local 
effectiveness  of  two  adjacent  measurement  regions,  28.5  >  N/p  >  27.0  and 
31.5  >  h/d  >30.  Secondary  flow  In  the  blade  passage  has  little  effect  on  the 
film  cooling  performance  at  distances  over  30  dlamters  from  the  end  wall. 

For  the  region  of  H/p  <  30,  the  variation  of  local  effectiveness  of  two 
adjacent  measurement  regions  Is  no  longer  periodic.  For  H/p  <  18,  local 
effectiveness  over  approximately  half  of  this  region  Is  under  0.1.  The  blade 
surface  Is  not  effectively  protected  by  the  film  cooling.  Much  lower 
effectiveness  Is  found  In  the  region,  H/p  <  1.5.  Two  reasons  are  believed  to 
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cmu  this  phenomenon.  First,  lower  mainstream  velocity  In  this  region 
results  In  a  higher  blowing  rote  end  coolant  blows  away  from  the  surface 
Instead  of  covering  the  surface.  Also,  the  hole  half  embedded  in  the  endwall 
Is  blocked  which  diminishes  the  coolant  flow  for  a  given  span. 

Due  to  limited  numbers  of  sampling  taps,  it  is  hard  to  determine  the 
downstream  location  at  which  the  local  effectiveness  just  reaches  zero. 

Because  the  locations  of  Iso-effectiveness  lines  are  calculated  by 
Interpolation,  the  Iso-effectiveness  lines  with  a  value  of  0.1  Involve  some 
degree  of  error. 

Figs.  4  and  5  show  the  spanwlse -eve rage  effectiveness  rt1w,  as  a  function 
of  the  dimensionless  distance  h/q-  Results  for  two  density  ratios  are 
presented.  Solid  symbols  on  the  right  axis  represent  <nTw  In  the  two- 
dimensional  flow  region,  at  H/p  •  69.  These  values  compare  well  with  results 
In  Ref.  [3].  The  data  at  each  downstream  position  (x),  have  a  similar  trend. 
As  H/0  decreases  the  average  effectiveness  at  first  deviates  only  slightly 
from  that  in  the  two-dimensional  flow  region.  At  a  certain  distance,  marked 
by  a  plus  symbol  In  Figs.  4  and  5,  the  values  of  the  average  effectiveness 
start  to  decrease  approximately  linearly  to  zero. 

For  each  specified  downstream  location,  the  results  can  be  approximated 
by  two  straight  lines.  One  extends  from  the  value  In  the  two-dimensional  flow 
region  and  Is  parallel  to  the  horizontal  axis.  The  slant  line  fits  data  which 
decreases  linearly  with  H.  Both  lines  Intersect  at  the  location  shown  by  a 
plus  symbol,  which  Indicates  where  the  value  of  the  average  effectiveness 
starts  to  decrease  linearly.  The  star  symbols  represent  the  locations  on  the 
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horizontal  axis  Intersected  by  the  slant  straight  lines.  The  region  between 
the  location  of  the  star  symbol  and  H/0  *  0  has  no  film  cooling  protection. 
Note  that  these  plus  and  star  symbols  are  not  data  points. 

The  plus  and  star  symbols  can  be  mapped  on  the  convex  surface  of  a  blade. 
Figure  6  presents  the  ranges  of  the  average  effectiveness  In  the  near  emfcall 
region  on  the  convex  surface  of  a  blade.  Note  that  the  location  at  a  value  of 
X/D  *  0  Is  at  the  downstream  edge  of  the  Injection  holes  of  the  second  row. 

The  leading  edge  of  a  blade  Is  at  x/g  -  -14.5.  The  star  symbols  (In  fig.  4  or 
5)  construct  the  upper  boundary  of  the  region  C.  The  upper  boundary  of  the 
region  B  Is  formed  by  the  corresponding  plus  symbols.  The  dashed  line  shows 
the  extropolated  boundary  thickness,  499,  of  the  mainstream  on  the  endwall  at 
the  leading  edge  of  a  blade. 

Region  C  Is  the  unprotected  area  where  coolant  flow  Is  swept  away  by  the 
passage  vortex  from  the  adjacent  blade.  The  span  of  this  unprotected  area 
grows  as  the  location  moves  toward  the  trailing  edge  because  the  passage 
vortex  becomes  bigger  as  It  travels  downstream.  In  region  B,  the  values  of 
the  average  effectiveness  vanish  linearly  to  zero  as  measurements  are 
conducted  toward  the  endwall.  In  this  region,  the  coolant  flow  mixes  with  the 
passage  vortex  but  Is  not  totally  absent.  Above  the  region  B,  the  performance 
of  the  film  cooling  Is  slightly  affected  by  skewed  streamlines  and  Is  close  to 
the  two-dimensional  flow  value  at  H/g  s  30. 

The  flow  field  near  the  end  wall  significantly  affects  the  performance  of 
the  film  cooling  on  the  convex  surface  of  a  blade.  The  passage  vortex  Is  so 
significant  that  a  large  area  of  a  blade  surface  could  be  exposed  to  hot  gas 
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without  coolant  protection.  It  may  be  helpful  to  add  coolant  to  the  end  wall 
to  protect  the  lower  convex  surface  of  the  blade.  It  Is  however  uncertain 
that  the  coolant  so  Injected  would  be  carried  up  by  the  passage  vortex  to  the 
blade  surface. 

The  effect  of  some  key  parameters  on  the  size  of  the  unprotected  region 
were  studied.  Fig.  6  shows  the  effect  of  the  density  ratio  on  boundaries 
between  the  regions.  At  a  blowing  rate  of  unity,  the  unprotected  area  Is 
smaller  with  the  lower  density  ratio  but  the  change  Is  slight.  Figs.  7  and  8 
show  the  effect  of  the  blowing  rate  on  the  three  regions.  Results  for  the  two 
density  ratios  studied  are  presented.  Higher  blowing  rate  might  be  expected 
to  decrease  the  size  of  the  unprotected  area.  However,  the  upper  boundary  of 
region  B  Is  not  affected  by  the  change  of  the  blowing  rate  and  there  Is  very 
little  affect  on  the  lower  boundary  of  B  (upper  boundary  of  region  C)  as  well. 
Figure  9  shows  the  effect  of  the  number  of  rows  of  Injection  holes  on  the 
regions.  Results  of  the  one-row  configuration  are  from  ref.  [4],  There  Is 
little  difference  between  the  two  sets  of  results  In  terms  of  the  size  of  the 
regions.  Outside  region  C,  the  values  of  the  average  effectiveness  of  the 
two-row  configuration  are  of  course  higher  than  those  for  one-row  Injection, 
but  the  locations  of  the  boundaries  are  not  very  different  for  the  different 
geometries. 

Concave  surface. 

For  M  ■  1.5  and  R  »  2.0,  Iso-effectiveness  lines  on  the  concave  surface 
are  shown  In  Fig.  10.  Measurement  regions  start  at  x/q  ■  1.10  and  end  at 
X/D  ■  38.86.  Because  streamlines  near  the  concave  surface  are  bent  toward  the 
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endwall,  periodicity  between  two  adjacent  measurement  regions  Is  not  present 
close  to  the  endwall.  The  Influence  of  the  endwall  on  local  effectiveness  can 
be  neglected  for  H/D  >  30.  Values  of  the  average  effectiveness  are  shown  In 
Figs.  11  and  12;  results  for  two  density  ratios  are  presented.  Values  of  the 
average  effectiveness  are  almost  Independent  of  H  except  very  close  to  the 
endwall.  The  presence  of  the  endwall  does  not  greatly  affect  the  performance 
of  the  film  cooling  on  the  concave  surface. 

Fig.  13  shows  a  comparison  of  the  average  effectiveness  between  one-row 
and  two-row  configurations  on  the  concave  surface.  Although  the  blowing  rate 
has  a  10%  difference,  a  comparison  Is  still  acceptable.  As  expected, 
effectiveness  with  two-row  Injection  Is  higher  than  with  one-row  Injection. 

In  the  region  H/p  <  3,  one  might  expect  that  the  values  of  the  average 
effectiveness  should  be  higher  than  those  of  the  two-dimensional  flow  region 
for  two  reasons.  The  first  Is  that  coolant  Is  washed  toward  the  endwall.  The 
other  Is  that  the  values  of  average  effectiveness  on  the  concave  surface 
generally  Increase  with  a  higher  blowing  rate  excluding  the  region  close  to 
the  Injection  hole.  From  Fig.  13,  this  Is  true  for  one -row  configuration  but 
not  for  two-row  configuration.  Due  to  the  blockage  of  the  lowest  hole  of  the 
second  row,  the  region  of  H/p  <  3  Is  not  fully  covered  by  two  rows  of 
Injection  holes. 

Conclusion 

The  presence  of  the  endWall  has  a  significant  effect  on  the  performance 
of  the  film  cooling  on  the  convex  surface  of  a  blade.  A  triangular  region, 
where  the  coolant  Is  swept  from  the  convex  surface  by  a  passage  vortex  from 


the  adjacent  blade,  has  no  film  cooling  protection  at  all.  Moving  away  from 
this  unprotected  area,  the  values  of  the  average  effectiveness  Increase 
linearly  from  zero  to  a  value  close  to  that  in  the  two-dimensional  mid-span 
region.  The  boundaries  of  these  regions  are  only  slightly  changed  by  varying 
blowing  rate,  density  ratio,  and  number  of  rows  of  injection  hole.  On  the 
concave  surface,  the  performance  of  the  film  cooling  is  not  significantly 
altered  by  the  presence  of  the  endwall  but  the  periodicity  of  the  local 
effectiveness  is  changed  due  to  skewing  of  the  streamlines  toward  the  endwall. 

The  film  cooling  effectiveness  is  closely  related  to  the  mainstream  flow 
field.  The  Incoming  flow  condition  and  the  geometry  of  the  turbine  blade  can 
play  Important  roles  in  affecting  the  dimensions  of  the  different  regions  of 
effectiveness  variation  on  the  convex  surface.  To  protect  the  region  where 
coolant  Is  washed  away.  It  may  be  necessary  to  further  cool  the  end  wall 
Itself. 
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mainstream  velocity  at  the  first  Injection  hole  location 
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density  of  foreign  gas 
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Fig.  13  Comparison  of  average  effectiveness  at  R  ■  0.96  between  one-row 
Injection  and  two -row  Injection 


